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INTRODUCTION 


For  solving  the  problem  to  be  treated  within  the  project,  it  was  advanced  the  mo¬ 
ment  method  for  solving  the  Boltzmann  equation.  Concerning  the  moment  method, 
the  distribution  function  (DF)  is  expanded  in  terms  of  the  Hermit  polynomials  Hj. 
Our  new  fundamental  result  obtained  several  years  ago  forms  the  basis  of  the  project, 
i.e.,  there  are  a  set  of  the  relationships  between  the  matrix  elements  (ME)  of  the 
collision  integral  [1],  [2], 

A  matrix  element  K3n  ]2  takes  a  sense  of  a  scalar  product  of  Hj  and  the  collision 
integral  of  Hj1  and  Hj2.  Namely,  the  lack  in  the  MEs  of  large  indices  retards  any 
development  of  the  moment  method  for  solving  the  Boltzmann  equation.  It  was 
shown  that  the  relationships  between  the  MEs  can  be  used  to  build  up  the  recurrence 
procedures  to  obtain  all  the  MEs  if  the  relatively  small  number  of  the  linear  MEs 
was  known.  This  is  the  very  value  of  the  new  result  to  be  applied  in  the  gas  kinetic 
theory.  It  should  be  noted  that,  in  the  past,  the  non-linear  MEs  were  known  in 
a  scarce  quantity  and  only  with  the  very  small  indices  and  the  formulas  for  the 
linear  MEs  (these  are  easily  expressed  via  the  integral  brackets)  were  known  only 
for  the  small  indices,  being  obtained  in  a  rather  general  form.  The  recurrence 
procedures  were  developed  and  built  up  to  calculate  the  non-  linear  MEs.  Within 
this  calculation,  the  linear  MEs  of  the  large  indices  were  built  up  to  be  its  integral 
part.  Now,  with  the  MEs  to  be  known,  it  should  be  maintained  whether  or  not 
the  moment  method  will  be  used  to  build  up  the  DF.  In  a  course  of  the  project 
work,  the  answer  ”Yes”  on  this  question  is  obtained.  Within  a  scope  of  this  project, 
the  linear  problems  were  treated  as  a  behaviour  of  a  small  ion  admixture  on  the 
atom  background  was  under  investigation.  Hence,  for  solving  these  problems,  it 
was  sufficient  the  linear  MEs  to  be  known.  However,  an  application  of  the  moment 
method  is  not  determined  by  the  linearity  or  non-linearity  of  the  collision  integral 
but  a  DF  form,  i.e.,  its  deviation  from  equilibrium.  The  strong  electric  held  is  one  of 
the  most  serious  sources  of  non-equilibrium  (in  the  case  of  the  charged  particles).  It 
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was  advanced  the  non-stationary  approach  which  was  successfully  used  for  solving  all 
the  problems  under  investigation  in  this  project.  The  essence  of  this  approach  is  in 
solution  of  a  problem  of  ion  evolution  after  a  sharp  switch-  on  of  the  electric  held.  As 
a  result,  not  only  the  stationary  state  is  found  but  the  details  of  the  non-stationary 
process  were  revealed.  The  developed  method  can  be  used  in  future  when  studying 
the  ion  behaviour  in  an  oscillated  electric  held.  The  problems  were  solved  for  the 
various  models  of  interaction  as  in  a  pure  electric  held  as  under  crossed  electric  and 
magnetic  helds.  In  the  latter  case,  a  series  of  interesting  phenomena  was  revealed 
on  the  non-stationary  stage  of  the  process.  All  calculations  were  carried  out  for 
the  case  of  equal  masses  of  the  ions  and  atoms.  Along  with  the  ion  mobility,  ion 
diffusion  was  studied  at  the  electric  held.  Here,  we  studied  successfully  two-stage 
non-stationary  approach:  the  DF  being  established  at  the  given  electric  held,  the 
non-stationary  process  was  studied  after  a  switch-on  of  a  density  gradient  which 
can  be  directed  as  along  the  held  (longitudinal  diffusion)  as  transversely  to  the  held 
(transverse  diffusion).  Treating  the  problems  inherent  for  the  project,  the  limits 
of  the  moment  method  application  were  found.  This  method  turns  out  to  be  used 
in  a  domain  of  the  moderate  and  not  so  strong  electric  helds.  In  a  region  of  the 
very  strong  electric  helds,  the  DF  moments  obtain  the  very  large  values  and  DF 
restoration  via  its  moments  can  not  be  accomplished.  All  that  is  due  to  the  Grad 
criterion  violation  for  application  of  the  moment  method  at  the  stationary  state.  We 
advanced  a  series  of  the  approaches  for  solving  this  problem,  i.e.,  the  DF  expansion 
in  terms  of  the  spherical  harmonics.  Here,  any  knowledge  on  the  MEs  with  large 
indices  gives  a  true  avenue  to  build  up  the  kernels  of  the  obtained  system  of  the 
integro-differential  equations. 

Our  work  on  the  project  was  proceeded  during  11  Quarters,  the  Quarter  12  being 
concentrated  on  this  Final  Report.  Below  is  the  list  of  the  problems  considered 
during  our  studies  on  the  project,  i.e.,  the  study  content  given  quarterly: 

Quarter  1.  Transfer  of  the  charged  particles.  Brief  literature  review.  A  choice 
of  the  way  of  the  study.  Matrix  elements  of  the  collision  integral.  New  results.  The 
Table  of  the  matrix  elements. 

Quarter  2.  A  system  of  the  moment  equations  in  different  bases. 

Quarter  3.  The  standard  moment  method  for  calculation  of  the  distribution 
function  and  ion  mobility.  Non-stationary  processes. 

Quarter  4.  Study  on  the  standard  method  convergence  for  the  ion  distribution 
function  calculation.  CEM-model. 

Quarter  5.  Numerical  solution  of  non-stationary  moment  system  for  various 
models  of  interaction. 

Quarter  6.  Comparison  of  the  ion  distribution  function  under  sharp  switch-on 
of  the  electric  held.  Details  of  the  moment  equation  for  the  charged  particles  under 
crossed  electric  and  magnetic  helds. 

Quarter  7.  Non-stationary  processes  for  the  charged  particles  under  mutually 
perpendicular  electric  and  magnetic  helds. 

Quarter  8.  Study  on  the  physical  moments  of  the  charged  particles  under 
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mutually  perpendicular  electric  and  magnetic  fields  for  various  models  of  interaction. 

Quarter  9.  Evolution  of  ion  distribution  function  under  mutually  perpendicular 
electric  and  magnetic  fields  at  a  sharp  switch-on  of  a  weak  electric  field. 

Quarter  10.  Evolution  of  ion  distribution  function  under  mutually  perpendic¬ 
ular  electric  and  magnetic  fields  at  a  sharp  switch-on  of  a  moderate  electric  held  for 
various  models  of  interaction. 

Quarter  11.  Study  on  ion  diffusion  under  the  electric  held. 

From  this  content,  one  can  see  that  every  point  which  follows  is  an  almost  nec¬ 
essary  continuation  from  the  preceding  ones;  it  is  difficult  to  change  this  sequence. 
At  the  same  time,  every  Quarterly  report  is  a  finite  solution  of  one  or  another  prob¬ 
lem.  The  Final  Report  is  to  be  logically  built  up  conserving  such  a  sequence  when 
formulating  in  brief  the  general  results  obtaining  in  the  quarters  from  1st  to  11th. 

In  the  Report  Appendix,  the  additional  results  obtained  for  the  Customer  request 
are  presented.  Here,  it  is  given  the  DF  presentation  as  well  as  that  of  ion  mobil¬ 
ity  under  the  electric  held  being  affected  simultaneously  by  the  resonance  charge 
exchange  and  isotropic  scattering. 

1  First  quarter 

First  quarterly  report  [3]  consists  of  three  parts. 

The  hrst  part  presents  a  brief  literature  review. 

There  are  a  lot  of  studies  on  mobility  of  the  charged  particles.  The  general  results 
of  the  investigations  accomplished  before  the  70s  of  the  last  century  can  be  found 
in  [4], [5], 

Last  decades,  the  most  advanced  studies  concerning  mobility  of  the  charged 
particles  have  been  performed  in  Australia,  K.  F.  Ness  and  R.  E.  Robson  being 
the  principal  authors  of  [6,  7,  8,  9].  In  Australia,  these  studies  were  undertaken 
on  Kumar  initiative  [10,  11,  12,  13]  and  continued  as  a  development  of  his  very 
interesting  article  [14], 

The  general  result  of  [14]  is  in  obtaining  the  direct  formulas  for  calculation  of 
the  matrix  elements  of  the  collision  integral.  In  a  general  case,  these  formulas  are 
very  cumbersome  and  useless  for  calculation  of  the  matrix  elements  with  large  indices 
because  of  the  multiple  sums  in  the  Talmi  coefficients.  However,  they  can  be  applied 
in  the  linear  ME  calculation. 

Most  of  the  general  results  were  obtained  for  the  electrons.  In  this  case,  the 
computation  problems  of  matrix  elements  are  easier.  From  this  literature  review,  a 
series  of  conclusions  was  deduced,  the  general  ones  being  as  follows 

1.  The  moment  method  is  the  key  when  solving  the  problem  which  we  are 
confronted  unless  the  unique  one. 

2.  The  main  problem  in  the  moment  method  is  calculation  of  matrix  elements 
of  collision  integral  with  large  indices. 
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3.  When  we  switch  to  the  weak  electric  field  from  the  strong  one  it’s  appro¬ 
priate  to  use  the  modified  moment  method.  In  this  method,  the  temperature  of 
weight  maxwcllian,  about  which  the  DF  expansion  is  undertaken,  differs  from  the 
temperature  of  atoms. 

4.  Potentialities  of  the  standard  moment  method  are  far  from  being  exhausted 
and  a  domain  of  its  application  has  not  been  specified,  the  method  of  computation 
of  any  arbitrary  matrix  element  being  lacking.  In  the  second  part,  the  new  results 

of  computation  of  matrix  elements  are  presented. 

The  moment  method  of  solving  Boltzmann  equation  is  based  on  the  polynomial 
expansion.  We  consider  the  polynomial  DF  expansion  in  terms  of  the  spherical 
Hermite  polynomials  [1],  [16]. 


/(v,r,t)  =  M(a,c)^Cj(r,t)Hj(c).  (1) 

3 

Hj  =  YL(Q,  v)clsl+1/2(c2),  i  =  0,1;  c  =  ^/a(v-u),  a  =  (2) 

<P)  =  pr(cos  0)  cos nrup  5^(0,  <p)  =  Pjm( cos  0)  sin rrup  (3) 

Here,  an  index  j  consists  of  four  indices  ( r,l,m,i ),  M  is  a  Maxwellian  distribution 
with  a  temperature  T  and  an  average  velocity  u,  and  Yim{6,  ip)  are  the  real  spherical 
harmonics,  .5'^,  ,2(c2)  are  the  Sonine  (Laguerre)  polynomials. 

The  Boltzmann  equation  is  replaced  by  a  system  of  moment  equations: 


D(Cj) 

Dt 


j,k 


(4) 


An  operator  D  /  Dt  in  (4)  corresponds  to  the  differential  (left)  side  of  the  Boltzmann 
equation.  When  there  is  not  any  external  force,  this  operator  was  studied  in  [16].  In 
the  presence  of  external  field  in  the  crossed  electric  and  magnetic  fields,  an  operator 
D/Dt  is  under  study  in  the  later  quarters  of  this  project. 

Poor  progress  in  the  non-linear  moment  method  is  due  to  the  large  difficulties 
when  calculating  the  MEs  which  are  determined  as  follows: 

/\'a.  =  /  H,I(MHV  MHk )  d3v/gz,  g,  =  f  MI  if  d3v,  (5) 

Here,  I  is  the  collision  integral,  M (T,  u)  is  a  Maxwellian,  and  g3  is  a  square  of  a 
norm  of  the  Hermite  polynomial  Hj. 

The  linear  matrix  elements  (LMEs)  of  the  types  K 0  and  Kq  J2  refer  to  the  linear 
collision  integral  and  could  be  of  the  first  type  Kju0  and  second  type  Kqj2.  Further, 
the  LMEs  are  denoted  as  A.  Thus,  e.g.,  in  an  axially  symmetric  case  we  have 

A  l1)  —  Kr’1  ■  A®  —  Kr'1 

J\,r i,l  ~  -rvn,Z,0,0>  1Vr,r2,l  ~  /Y0,0,r2,/  VU/ 
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When  calculating  the  kinetics  of  an  ion  admixture,  we  need  the  LME  of  the  first 
kind. 

In  the  ID  spatial  problems,  the  velocity  distribution  function  is  the  axially  sym¬ 
metric  one,  and  it’s  expansion  is  performed  in  terms  of  the  spherical  Hermite  poly¬ 
nomials  of  two  indices  (r,l):  Hlr  =  Sf+1/2Pi. 

In  [1],  it  was  shown  that  3 D  MEs  are  proportional  to  the  corresponding  axially 
symmetric  MEs  and  the  coefficients  of  proportionality  are  easily  expressed  via  the 
Klebsh-Gordan  coefficients. 

For  convergence  of  (1),  the  distribution  function  should  satisfy  the  Grad  criterion 

poo 

/  /2exp(c2)d3u  <  oo.  (7) 

Jo 

The  limits  concerning  the  Grad  criterion  arise  when  an  expansion  is  Carried  out 
in  terms  of  the  Sonine  polynomials  but  not  on  the  stage  of  the  expansion  in  terms  of 
the  spherical  harmonics  because  the  Sonine  polynomials  do  be  the  orthogonal  ones 
with  the  Maxwellian  weight. 

As  a  base  of  new  relationships  between  the  MEs,  we  advance  the  collision  inte¬ 
gral  invariance  relative  on  a  set  of  the  basic  functions,  in  terms  of  which  the  DF 
is  expanded.  The  bases  distinct  in  the  parameters  of  the  weight  Maxwellian  are 
considered.  The  coefficients  of  the  DF  expansion  are  interrelated  via  a  transfer 
matrix 


C1  =  Ebm(w'i,W'o)C2.  (8) 

k= 0 

In  the  axially  symmetric  case,  W  incorporates  two  parameters  T  and  u.  Matrix  D 
was  built  up  in  [1]  with  the  use  of  a  —  u  presentation  of  the  Boltzmann  equation. 

An  important  theorem  results  from  the  invariance  principle  of  the  collision  inte¬ 
gral  relative  on  a  choice  of  the  basic  functions:  the  MEs  of  the  collision  integral  in 
two  different  bases  are  expressed  via  each  other  as  follows 


KAW i)  =  E  Di'AW 1>  Wo)  53  Klj(W0)Dk  kr(W0,  WjDjAWo,  WJ.  (9) 

i  k,j 

This  expression  being  differentiating  with  respect  to  T\  or  u \  and  setting  W±  = 
Wo  ( Ui  =  Uq,  T\  =  T0),  the  following  relationships  between  the  matrix  elements 
arise. 

”  Temperature”  recurrence  relationships 


M'( 


a  \ 
a,  b) 


r,l 

ri,h,r2,h 


RK{ 


a  \r,l 

a,b)ri,h,r2,h 


+  rK(i  b) 


r—l,l 

ri,h,r2,l2 


-(r1  +  l)W(“fe) 


r,l 

ri+l,h>r2,h 


(r2  +  l)K(l  b) 


r,l 

n,Ji,r2+i,/2  ’ 
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(10) 


Here,  R  —  r i  +  r2  —  r  +  (h  +  l2  -  0/2- 

In  the  case  of  the  exponent  potentials  (V  ~  -77),  the  MEs  are  proportional  to  TM, 
where,  e.g.,  fi  =  0  for  the  Maxwellian  molecules,  fi  =  0.5  for  hard  spheres,  /i  =  —1.5 
in  the  case  of  Coulomb  interaction. 

”  Velocity”  recurrence  relationships 


0 


M  -  1  maa,b) 


r,l—  1 
ri,h,r2,h 


+  7(r 


i,(  +  i)*-(y 


r— 1,/+1 
ri,h,r2,l2 


-mK(aa,b 


r,l 

ri,li+l,r2,h 


r,l 

ri+l,h  —  l?r2^2 


(ii) 


-yj(mh/ma)  (j3{l2)K{aa  h) 


r,l 

ri,h,r2,l2+l 


7  Mmit. 


1*7 


/?(/) 


/+  1  7(r  /)  =  (r  +  ^ 
2/  +  1  ’  71  ’  ’  21  +  1 


Here,  we  consider  a  mixture  of  particles  of  two  kinds,  a  and  b,  with  the  masses  ma 
and  mb-,  correspondingly.  The  equations  are  written  for  the  ME  of  a-kind,  scattered 
by  the  particles  of  6-kind.  The  recurrent  equations  (10)  and  (11)  hold  not  only  for 
the  matrix  elements  of  the  full  collision  integral  but  separately  for  the  MEs  of  the 
incoming  and  outgoing  terms. 


For  the  ordinary  (not  oriented)  particles,  the  linear  Boltzmann  operator  meets 
the  Hecke  theorem  [18],  which,  in  terms  of  the  matrix  elements,  says:  K^li00  =  0 

and  Ko!o,r2,i2  =  0,  if  and  l2  ^  l. 

Using  a  recurrence  formula  (10),  (11),  it’s  easy  to  prove  the  generalized  Hecke 
theorem. 


* 0,  if  \h-i2\<i<h+k  (12) 

and  a  parity  of  /  coincides  with  a  parity  of  /1  +  l2.  The  recurrence  procedures  were 
developed  and  it  was  shown  that  all  the  non-linear  (and  linear)  MEs  could  be 
found  if  the  linear  isotropic  (1=0)  MEs  of  the  same  type  are  known.  The 

simple  formulas  [1],  [2]  were  obtained  for  the  linear  matrix  elements  of  the  second 
kind. 

A  program  package  which  allows  to  compute  all  MEs  including  the  MEs  with 
the  very  large  indices  for  arbitrary  masses  of  colliding  particles.  Using  these  pro¬ 
grams,  the  properties  of  the  MEs  were  studied  and  some  interesting  features  in  their 
behaviour  were  found. 

New  obtained  results  open  new  potentialities  of  the  moment  method  for  compu¬ 
tation  of  the  DF  in  a  domain  of  the  high  energies. 

In  the  third  part  of  report,  the  Tables  of  the  MEs  for  different  models  of  inter¬ 
action  were  presented.  A  comparison  with  the  known  literature  data  were  made. 
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2  Second  quarter 


Second  quarter  was  dedicated  to  an  investigation  of  a  differential  operator  of  the 
left  side  of  the  Boltzmann  equation  [19]. 

When  solving  the  spatially  non-uniform  problems  by  the  moment  method,  a 
basic  Maxwellian  can  be  chosen  as  the  unique  one  in  all  the  spatial-temporal  points. 
Such  a  Maxwellian  is  name  as  the  global  one.  Other  choice  is  a  local  basis,  when, 
for  any  point,  a  velocity  of  the  weight  Maxwellian  coincides  with  a  mean  velocity 
of  a  gas  in  this  point,  and  a  temperature  is  determined  via  a  mean  thermal  energy 
of  the  gas.  Note,  that,  besides  of  the  global  and  local  bases,  the  other  ones  can  be 
chosen,  being  the  intermediate  ones  between  the  formers. 

When  building  up  the  left  side  of  the  moment  equations  in  the  global  basis,  the 
Boltzmann  equation  is  suitable  expressed  as  follows 

Af  =  /(/,  /) 

Here,  /(v,  r,  t)  -  is  the  DF,  /(/,  /)  -  is  the  collision  integral. 

2  r  _  df  ,  df  F  df 

Af  =  y:  +  v  ■  77 — I - -  yr"  > 

at  or  mo  ov 


(13) 


A  F  is  a  vector  of  the  external  force,  a  mo  is  a  particle  mass.  An  operator  A  is  a 
linear  differential  operator  of  the  left  side  of  the  Boltzmann  equation. 

In  the  method  of  the  polynomial  expansion  in  terms  of  the  spherical  Hcrmite 
polynomials,  the  DF  is  presented  as  (1).  The  expansion  coefficients  Cj  are  expressed 
via  the  DF 

Cj(r,t)  =  -J-  f  f(v)Hj(c)d3v.  (15) 

9j  J 

A  normalizing  factor  gj  is  as  in  [1] 


27t(1  +  5mo  )(l  +  m)\  r(r  +  l  +  3/2) 

(l  —  m)\(2l  + 1)  ’  rl  27r3/2r! 


(16) 


The  DF  is  substituted  of  its  equivalent  expression  via  an  infinity  vector  C  with 
the  components  (15).  Denote  the  left-side  operator  in  this  expression  via  P.  Then, 
we  have 

Q  =  PC.  (17) 

The  components  of  a  vector  Q  are  the  left  sides  of  the  corresponding  moment 
equations.  Rewrite  (17) 

Qj  =  (18) 

y 

The  elements  Pjy  are  expressed  via  the  MEs  of  an  operator  A  as  follows 


Pjy  =<j\A\j'>, 
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where 


<j\A\f  >=  —  /  Hj(c)A(M(a,c)Hr(c))  d3v. 
9i 


Now,  we  have 


or 


Qi  =  DCj 


3  -  ~JZr  t-L  : ,4/)  +«1/2  (— ,ic/ 

Dt  a1/ 2  V  .9,  ;  V  gj  , 


DC,  dCj 


Dt 


dt 


+  u 


<9r 


where 


+/  =  c 


d£ 

dr' 


Aj  = 

m0  d  c 


A  scalar  product  (Hj/gj,  Arf)  can  be  expressed  as 


—  >  ^r/1  =  -  /  ^i(c)  +  Cy^  +  c/+. 


■A 


9j 


'  dz 


or 


'ff,-  ,  A  3  9 


->4/  =E^f, 

9j  J  k=l  dxk 

G[k)  =  -  [  ckHjfd3v. 

9j 


After  the  cumbersome  operations,  the  expressions  which  follow  are  obtained 
Gj  "*  =  —  2^ j  _|_  gj  A  r  +  3/2)dm_iCrZ+lm_1  —  rdm_iC,r_1;+lm_1 
-(/  +  m  +  2) (l  +  m  +  1)  ((/  +  r  +  3/2 )Clrl+lm+1  —  rC* 

— lZ+lm+l)  ^ 

i^dm.—  liGirl—lm—l  ^r+lZ— lm—  l) 


2(2/-  1) 

—  (/  —  m)(l  —  m  —  1)  (C, 


G 


rl— lm+1  '“/r+l/— lm+1  y  J  ? 


(19) 

(20) 

(21) 

(22) 

(23) 

(24) 

(25) 

(26) 


(27) 


Gf  = 


■1) 


i+1 


(/  +  r  +  3/2 )  c/m_  i  Gw+lm_  x  rdm_iGr_1;+lm_1 


2(2/ +  3) 

+  (/  +  m  +  2)(/  +  m  +  1)  ((/  +  r  +  3/2)G/;_/lm+1  —  rC/++lm+1 
f_iV+1  / 

V  /  (  ]  (/"!  l—i  _  y^fl—  i  \ 

»i-l  v'-'rZ— lm— 1  '-yr+lZ— lm—  1/ 


2(2/-  1) 

+  (/  -  m)(/  -  m  -  1)  (C3t 


—i  _  ri\—i 

rl— lm+1  ^r+lZ— lm+1  /  I  ? 


(28) 


G?  =  -  \  ™ ' 3  1  ((l  +  r  +  3/2 )C'rl+lm  -  rC‘_ll+lm) 
l  —  m  /  ■  •  \ 


2^  ^  Ira  wr+lZ— lm J  ' 


Here,  we  use  a  denomination 


dm  1  V  V, 


For  the  terms  related  with  the  external  force,  we  have 

(—,Af)  =-f  H,(  c)  t,  — |Vv  (31) 

\9j  J  9j  '  mo  dck 

Via  integration  by  parts,  taking  the  very  fact  that  the  DF  turns  out  to  be  zero 
at  each  infinity  limit,  (31)  can  be  expressed  as  follows 


9j  J 


Lf  -  f  fJ^—HjCpv. 

9j  J  ock 

The  following  expressions  for  dE-dE  are  obtained 


T  (1)  _ 
L3  ~ 


21  +  3 


21  -  1 


(/  -  m){l  -m-  1  )Clrl_lm+1  - 


(I  +  m  +  1)(/  +  m  +  2 )Clr_ll+lm+1  -  Clr_ll+lm_1 ) , 


'  (2)  _  /  1  Yi+1  ' 


Lj  =  (-1  )* 


21-1 


((*  -  m)(l  ~m~  l)Cll- -lm+1  +  C'w-lm-l) 


21  +  3 


(l  +  m+  1  (/  +  m  +  2  )Cr_[l+lm+l  +  Hr_[l+lm_l  ] , 


L f  =  y  +  +  m  +  W-lI+lm  “  2[ri(*  -  m)Crt-lrn-  (36) 

Thus,  finally,  the  left  side  of  the  moment  system  with  the  force  terms  in  the 
global  basis  is  as  follows 


Qi  = 


do.,  d  a 


4  E  4  °jk)  +  «1/2  E  — L 

1/1/2  <  J  ^  rm  „ 


dr  a1/2 


S  mo 


The  left-side  operator  of  the  moment  system  can  be  expressed  as  follows 

D  1  „  d  a1/2 

P  =  -T v  +  ^(r't)  +  — (F'S)- 

Ilf  a1/2  dr  m0 
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Any  specific  form  of  the  MEs  of  operators  R  and  S  is  easily  found  from  the  formulas 
for  Gj  (27)-(29)  and  Lj  (34)-(36),  respectively. 

When  deducing  the  left  side  of  the  moment  system  in  the  local  basis,  a  relation 
of  the  expansion  coefficients  the  DF  Cj  in  the  different  bases  is  used  via  a  transfer 
matrix  D  (8). 

The  basis  is  featured  by  a  4-vector  W,  of  which  components  are  the  Maxwellian 
parameters:  u  is  a  mean  velocity  and  T  is  a  temperature. 

When  using  a  local  basis,  W  depends  on  a  coordinate  and  time.  When  deducing 
an  equation  in  a  point  (ro,£o)  in  the  local  basis,  we  take  in  mind  the  equations  in 
the  global  basis  with  4-vector  W,  which  coincides  with  W  of  the  local  basis  in  this 
point.  In  a  neighbourhood  of  a  point  ,  (r0,  £o)a  relation  between  the  vector  C  in  the 
global  frame  (Cal)  and  the  local  one  ( Cloc )  is  as  follows 

C  loc(Wloc)  =  D(Wloc ,  Wgl)Ca\Wgl).  (39) 

In  a  point  r0,£o;  we  have  Wioc  =  Wgi,  and  a  transfer  operator  from  the  global 
basis  to  the  local  one  D(Wioc,  Wgi)  turns  out  to  be  an  identity  operator.  A  linear 
operator  P  is  transformed  at  the  transfer  from  the  global  basis  to  the  local  one  as 
follows 

Ploc  =  D(Wloc ,  Wgl)Pal(D{Wloc,  Wgl ))"1  =  D(WloC)  Wgi)P9lD(Wgl,  Wloc).  (40) 

For  the  left  sides  of  the  system  of  the  moment  equations  we  have 

Qioc  =  piocQioc .  (41) 

An  operator  Pal  is  a  differential  operator.  So,  at  Wioc  =  Wgi,  one  can  replace,  in 
(40),  D(Wioc,  Wgi),  at  the  left,  by  the  identity  one  at  WXoc  =  Wgh  but  D(Wgh  WXoc ) 
under  a  sign  of  the  differential  operator  can  not  be  excluded.  Now,  from  (41),  using 
(40),  we  obtain 

Q'“  =  {P*b( WlU  W,o0) C'“) \w„,-w,r  (42) 

As  an  operator  Pal  incorporates  the  derivatives  with  respect  to  r  and  t,  and 
the  transfer  operator  depends  on  these  variables  via  Wioc,  now,  the  spatial-temporal 
derivatives  should  be  applied  to  a  product  D(Wgi,Wioc)Cloc,  and,  taking  in  mind 
(38),  where  a  Pal  is  determined,  we  have 

Q'“  =  P* C'“  +  {  ((A  +  -lj(R  ■  B(W,U  WuS)  !_„,.}  C'“.  (43) 

In  the  local  basis,  the  second  terms  arises  in  addition.  Denote  it  as  Q.  Take  attention 
to  a  lack,  here,  the  force  terms  as,  in  this  part  of  an  operator  Pa\  there  are  not  any 
spatial  and  temporal  derivatives. 
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After  a  series  of  transformations,  including  a  search  of  the  derivatives  of  a  matrix 
D  with  respect  to  u  and  T,  we  obtain 


Qf  =  “Trh0)  +  <*1/2  £ 


Dt 


k= 1 


Dt 


(44) 


In  the  axially  symmetric  case,  i  —  m  —  0,  and  there  is  only  unique  component 
of  a  vector  u,  denoted  by  w.  An  equation  (44)  takes  a  form 


Here, 


Qf  =  a 


£Z>)  +  al/2^L(3) 
Dt  J  Dt  J 


>dT  ^  _  ,(0)  <9w 


(3) 


^  i'  ^  i' 


(45) 


jT  =  (r + 


r  (3)  _ 

W  - 


(ioc)  _  Mice) 
rCV-l/> 

r(l  +  1)C<‘“|+1  - 


2/ +  3 

(0,3)  l  +  1 


-ia 


(loc) 


rl—1  ? 


\  A  r  (U,3)  _ 

2_^  njj,L/j' 


2/ +  3 


2/  -  1 

((/  +  r  +  3/2)LS;3i  -  rLj°jS+1 


,  l  /  r  (0,3)  r  (0,3)  A 


(46) 


We  verified  that  the  continuity,  motion  and  energy  equations  in  a  form  advanced 
in  [20]  are  results  from  several  first  moment  equations  in  the  local  basis. 


3  Third  quarter 

In  Quarter  3  [21],  the  research  ways  to  be  developed  during  several  next  quarters 
were  formulated  concerning  ion  behaviour  under  the  external  electric  field. 

The  project  as  a  whole  is  devoted  to  computation  of  the  transport  coefficients 
and  the  DF  of  an  ion  admixture  on  an  equilibrium  (Maxwellian)  gas  background  via 
the  moment  method  either  under  the  external  electric  field  or  crossed  electric  and 
magnetic  fields. 

In  those  problems,  the  collision  integral  is  linear,  and,  in  the  moment  method, 
it  corresponds  to  the  linear  MEs  of  the  first  kind.  Note,  that  the  linearity  of  the 
collision  integral  does  not  exclude  a  strong  ion  non-equilibrium.  Under  an  increase 
in  the  electric  field  strength,  a  deviation  from  the  Maxwellian  distribution  could  be 
very  large.  Hence,  the  general  problems  of  the  moment  method  convergence  can  be 
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studied  with  these  (linear)  problems  as  an  example.  When  working  on  this  project, 
we  consider  the  ions  of  which  mass  coincides  with  a  mass  of  atoms.  We  denote  the 
standard  moment  method  as  that  in  which  the  distribution  function  expansion  is 
carried  out  about  the  Maxwellian  with  a  temperature  of  atoms.  Along  with  that, 
an  expansion  can  be  performed  about  other  Maxwellian  of  which  parameters  do  not 
coincide  with  the  parameters  of  the  background  atoms  distribution;  such  approach 
is  denoted  as  the  modified  moment  method.  Under  moderate  fields,  the  standard 
moment  method  was  consider  by  Kihara  [22],  He  developed  an  iteration  method  to 
be  applied  with  the  minor  modifications  in  [4],  Any  knowledge  on  the  elements  of  the 
collision  integral  is  urgently  needed  to  succeed  in  the  moment  method  application. 
With  an  increase  in  the  electric  held  strength,  the  number  of  the  MEs  needed  for 
the  iteration  process  convergence.  In  [4],  the  formulas  for  the  ME  calculation  at  no 
so  large  indices  were  tabulated.  The  authors  emphasized  that  this  collection  of  the 
formulas  is  the  most  general  one,  incorporating  the  Ferziger,  Caper  results  [23]  as 
well  as  the  authors  themselves.  Overall,  there  are  the  formulas  for  calculation  of 
about  two  tens  of  the  MEs  in  these  Tables. 

We  compared  the  MEs  calculated  via  the  formulas  of  [4]  with  our  results  shown 
in  the  Tables  in  the  first  Quarterly  report.  It  turns  out  that  there  is  a  series  of 
distinctions  in  the  ME  definitions:  opposite  sign  in  the  collision  integral,  differ¬ 
ent  definition  of  the  expansion  coefficients  of  the  DF  and  different  choice  in  of  a 
measurement  unit  of  collision  frequency.  Further  comparison  results  in  an  ultimate 
coincidence  of  our  data  with  calculation  via  the  Tables  from  [4]  except  a  single  ME 
because  of,  eventually,  with  misprint  in  the  corresponding  formula  in  [4],  taking  in 
mind  all  the  distinctions  mentioned  above. 

To  perform  the  calculations,  a  transfer  to  the  dimensionless  variables  was  carried 
out.  As  a  velocity  unit  it  was  chosen  always  the  thermal  velocity  of  atoms. 


As  a  time  unit,  that  of  between  collisions,  r,  was  chosen.  If  a  scattering  cross  section 
~  1/gA,  where  g  is  a  relative  velocity  of  the  colliding  particles,  r  does  not  depend 
on  a  velocity;  a  constant  frequency  between  collisions  (a  time  between  collisions) 
is  inherent  for  the  Maxwellian  and  similar  to  the  Maxwellian  (pseudo-Maxwellian 
molecules;  BGK-model).  In  these  cases,  a  choice  in  time  unit  is  clear,  and,  when 
transferring  to  these  units,  for  dimensionless  electric  held,  we  obtain 

eEri’T  _  eEr 
m  V2  kTm 

However,  for  other  interaction  cross  sections,  the  collision  frequency  depends  on 
a  relative  velocity,  and  a  choice  in  time  unit  is  not  so  clear.  E.g.,  for  hard-sphere 
model,  the  scattering  cross  section  does  not  depend  on  a  velocity,  and  it  turns  out 
that  not  the  collision  frequency  is  constant  but  a  free  path  length  A. 

A  =  Nand\  (49) 
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Where  Na  is  the  concentration  of  atoms,  and  d  is  a  diameter  (in  general  case,  a 
sum  of  radii  of  colliding  particles). 

The  MEs  are  calculated,  also,  at  the  definite  normalization,  and  always  the 
unique  question  arises  what  to  be  included  in  a  definition  of  collision  frequency  and 
what  should  be  incorporated  in  the  ME.  In  [4],  the  collision  frequency  for  hard 
spheres  of  equal  masses  is  chosen  as  follows 

J-  =  a  (iWf 
tm  V  nm  ) 


When  calculating  the  dimensionless  MEs,  we  choice  the  collision  frequency  as  follows 


1 

te 


(51) 


In  these  very  units,  in  the  first  quarter,  the  ME  Tables  were  presented.  In  the  third 
quarter,  as  the  time  unit,  it  was  advanced  te,  when  transferring  to  the  dimensionless 
Boltzmann  equation,  and  such  normalization  was  named  the  standard  one.  In  these 
units,  for  the  dimensionless  electric  held,  we  have 

eEX 
V2kT 

Further,  beginning  in  the  sixth  quarter,  a  definition  of  the  standard  normalization 
was  changed;  conserving  the  velocity  unit  (47),  a  time  unit  was  changed.  For  the 
interaction  models  with  a  constant  free  path  length,  naturally,  A  should  be  chosen  as 
a  length  unit.  With  the  chosen  units  of  length  and  velocity,  r  is  determined  uniquely 
as  r  =  X/vt ■  Thus,  for  the  particles  with  the  velocity  independent  scattering  cross 
section,  the  standard  normalization  was  named  that,  for  which 

i=a(— y/2 

r  \  m  J 


(52) 


When  transferring  to  the  final  standard  normalization,  we  need  to  denominate  again 
all  the  dimensionless  MEs,  namely,  to  divide  by  y/2  for  the  models  with  A  =  const. 
For  the  final  standard  normalization  for  e  we  have 


_  eEX 
6  ~  2kT' 

As  it  was  noted  earlier,  all  calculations  beginning  from  the  sixth  quarters  are 
presented  in  the  final  standard  normalization.  However,  in  the  fifth  quarter,  for  the 
models  with  A  =  const,  in  the  results  involving  the  standard  normalization,  one 
should  denominate  all  the  values  of  e  and  t,  namely,  e  is  to  be  divide  by  y/2,  and  a 
time  is  to  be  multiplied  by  y/2,  referring  to  the  Figure  captions  as  well  as  in  text. 
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Go  to  the  dimensionless  variables.  As  a  velocity  unit,  we  choice  a  thermal  velocity 
vt  and  denote  a  dimensionless  velocity  c 


c  = 


V 

vt  ’ 


(55) 


and,  for  the  dimensionless  DF  and  time,  we  conserve  the  previous  denominations: 
/(c,f),  and  t.  In  the  dimensionless  form,  the  Boltzmann  equation  takes  a  form 


d£  df_ 

dt  e  dcz 


(56) 


where  J(/,  M )  is  the  dimensionless  collision  integral.  In  an  equation  (56),  the  z-axis 
direction  coincides  with  the  electric  held  direction.  The  corresponding  system  of  the 
dimentionless  moment  equations  takes  a  form 


dCrj 

dt 


+  6 


2 

— - r 

2/ +  3 


(/  +  l)C,y_i;_|_l 


21 

21  -  1 


Cr,l- 1)  = 


—  ^2 

7*1 


y  7*1,1 


(57) 


Here,  if  the  terms  with  a  temporal  derivative  are  omitted  and  the  distinctions 
in  the  definitions  of  the  collision  frequency  as  well  as  normalizations  are  taken  into 
consideration,  the  system  coincides  as  whole  with  the  stationary  system  in  [4],  for 
which  solution,  Kihara  [22]  advanced  the  iteration  method. 

As  a  rule,  when  solving  the  kinetic  equations  for  the  ions  at  the  external  held,  the 
stationary  state  of  a  system  is  looked  for,  i.e.,  the  stationary  values  of  mobility  and 
the  stationary  DF.  Namely,  in  the  third  quarter,  we  advanced  the  non-stationary 
approach  to  find  the  stationary  solutions.  For  brevity,  we  call  such  approach  the 
non-stationary  moment  method.  In  this  method,  the  full  system  of  the  moment 
equation  is  used,  including  the  terms  with  a  temporal  derivative,  and  it  is  supposed 
that,  till  7  =  0,  the  system  is  in  a  total  equilibrium,  i.e.,  the  ions  has  the  Maxwellian 
distribution  coinciding  with  the  atom  one.  At  t  —  0,  the  electric  held  is  instantly 
switched  on.  By  solving  the  non-stationary  system  of  the  moment  equation,  the 
transient  process  is  studied,  beginning  with  the  electric  held  switch-on  and  termi¬ 
nating  when  obtaining  the  stationary  state.  As  a  result,  we  find  the  ion  mobility 
and  DF  at  the  stationary  state  with  no  use  of  the  complicated  iteration  procedure. 
Besides,  the  ion  behaviour  in  a  course  of  the  transient  process  is  of  interest  itself. 

The  non-stationary  approach  was  used  successfully  during  all  the  work  on  this 
project  including  the  study  on  the  effect  of  the  crossed  electric  and  magnetic  fields 
on  the  ions  as  well  as  the  ion  diffusion  under  the  electric  held. 

Clearly,  to  solve  via  the  moment  method  the  problems  formulated  above,  the 
number  the  involved  moment  equations  should  be  essentially  larger  as  compared 
with  the  stationary  mobility  calculation. 

Development  of  an  algorithm  to  solve  the  non-stationary  problem  as  a  perspective 
gives  an  opportunity  to  go  to  an  important  applied  problem  concerning  the  ion 
behaviour  under  alternate  electric  held. 
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When  solving  the  non-stationary  moment  system,  the  4th-order  Rnnge-Kutta 
method  is  used. 

When  studying  the  moment  method  convergence,  it  is  of  importance  to  have  an 
analytical  solution  of  the  problem  under  consideration,  at  least,  in  several  particular 
cases.  Such  an  interaction  cross  section  was  obtained  with  which  an  analytical 
solution  not  only  of  the  stationary  problem  but  the  non-stationary  process  could  be 
built  up. 

Under  consideration  is  a  model  in  which  a  cross  section  is  inversely  proportional 
to  the  relative  velocity,  and  the  angulat  part  takes  a  form  8(9  —  n).  This  stationary 
model  of  the  charge  exchange  with  the  Maxwellian  dependence  of  the  cross  section 
on  the  relative  velocity  we  call  the  CEM-model.  As  the  collision  frequency  in  this 
model  does  not  depend  on  the  velocity,  it  is  eliminated  from  the  collision  integral; 
the  outgoing  term  of  the  collision  integral  equals  — /(v)/t)  where  /(v)  is  the  ion  DF. 
In  it’s  turn,  (5-shaped  angular  dependence  results  in  proportionality  of  the  incoming 
term  to  the  atom  DF.  With  a  small  ion  admixture  on  the  atom  background,  the 
atom  DF  is  in  equilibrium,  i.e. ,  it  takes  a  form  of  the  Maxwellian  distribution.  As 
a  result,  for  GEM-  model,  a  collision  operator  takes  a  form 

i  =  M(C)-/(C).  (58) 

T 

In  it’s  form,  it  coincides  with  the  linear  BGK-model  [24],  which  is  frequently  used  to 
obtain  the  analytical  and  numerical  solutions  in  a  gas  of  unique  sort.  If  the  BGK- 
model  is  made  with  aim  to  simplify  the  collision  integral  and  has  not  the  physical 
ground,  on  contrary,  the  CEM-model  is  an  implication  of  a  choice  of  a  well  defined 
scattering  cross  section.  Simultaneously,  it  has  all  the  gains  of  BGK-model  and  can 
be  used  to  build  up  the  analytical  solutions.  Present  with  no  derivation  the  general 
formulas  describing  the  analytical  solutions.  At  the  chosen  measurement  units,  a 
drift  velocity  of  ions  equals 

vd  =  C( 0,  l)/2, 

where 

C(0, 1)  =  2e(l  —  exp(-t)).  (59) 

Thus,  in  the  chosen  measurement  units,  the  mobility  goes  to  unity  at  t  — >  oo. 

The  stationary  solution  for  the  DF  takes  a  form: 

f(cP,cz )  =  ^-M(cp)M(cz)^~  exp  (cz  -  (l+erf(cz  -  .  (60) 

The  non-stationary  solution  for  the  DF  at  the  axis  of  symmetry  takes  a  form: 

f(t,cz)  =  ^M(c^)e(c2"1/(2e))2i  (erf(et  -  (cz  -  ^-))  + 

+erf(cz  -  +  e~tM(cz  -  et).  (61) 
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To  obtain  these  solutions,  the  method  of  characteristics  was  used.  Note,  that,  for 
CEM-model,  the  solution  can  be  built  up,  also,  for  other  non-  stationary  processes, 
e.g.,  to  describe  the  ion  behaviour  under  periodical  external  electric  held.  However, 
a  solution  of  such  a  problem  is  far  from  the  scope  of  this  report. 

As  an  example,  Fig.  1  shows  the  calculation  results  via  the  analytical  formulas 
of  the  DF  at  e  =  2.  Figure  la  is  the  DF  dependence  on  Cz  at  the  axis  of  symmetry 
(cp  =  0)  at  the  different  times.  Figure  lb  is  an  analogous  dependence  referred  to  the 
stationary  DF.  Here,  it  is  introduced  a  denomination  fr( 0,  cz)  =  /( 0,  cz)/ fstac( 0,  cz ), 
where  fstac( 0,  cz )  is  the  DF  at  the  axis  of  symmetry  at  the  stationary  state.  At 
t  — >  oo  fr  — >  1.  In  a  velocity  space,  along  the  electric  held,  one  can  see  the  motion 
of  a  front  of  which  behind  the  stationary  DF  arises. 


Figure  1:  Distribution  function  f(cz,  t )  at  the  axis  of  symmetry  at  the  different  times 
after  a  switch-on  of  the  the  electric  held  strength  of  e  =  2  (a)  and  the  distribution 
function  referred  to  the  stationary  DF  fr(cZlt)  =  f(cz,  t)/f(cz,  oo)(b).  The  times 
t  =  0  correspond  to  curve  (1),  0.5  (2),  1  (3),  2  (4),  3  (5),  4  (6),  5  (7),  10  (8). 

As  it  can  be  seen  in  Figure,  relaxation  in  a  domain  of  the  large  positive  values 
of  cz  proceeds  exclusively  slow,  and  the  stationary  solution  yielded  in  this  velocity 
domain  developed  with  a  large  delay. 

When  studying  the  moment  method  convergence,  the  Grad  criterion  feasibility  is 
of  crucial  value,  for  which  aim  a  sufficiently  rapid  decrease  in  the  DF  at  the  infinity 
is  needed.  It  is  shown  that  the  stationary  solution  for  the  interaction  model  under 
consideration  does  not  satisfy  the  Grad  criterion.  However,  at  any  finite  time  t  in 
a  course  of  the  non-stationary  process,  the  DF  decreases  in  a  large  velocity  domain 
in  a  such  rapid  manner  that  the  Grad  criterion  holds  and  there  are  not  any  general 
limits  to  the  moment  method  convergence.  As  a  result,  in  a  course  of  the  non- 
stationary  process,  a  large  part  of  the  stationary  DF  can  be  built  up.  Besides  of  the 
analytical  solution  for  the  DF  in  the  case  of  CEM-model,  we  obtain  successfully  the 
analytical  solution  for  the  moment  of  the  DF. 

Crtl{t)  =  BrjeN  SN(t)t  N  —  l  +  2  r.  (62) 
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Here,  Br  l  is  a  coefficient  independent  on  either  time  or  e.  It  equals 


(2;  +  4r)!!(2i  +  l) 

'  2r(2r  +  21  +  1)!! 


A  function  5A(f)  takes  a  form 


(63) 


N- 1  ik 

SN(t )  =  1  -  exp(-t)  ^2  — ,  (64) 

k= 0  K~ 

and  it  is  the  same  for  total  layer  N  =  l  +  2 r.  This  function  turns  out  to  be  zero  at 
t  —  0  and  1  at  t  — >  oo. 


4  Fourth  quarter 

Quarter  4  deals  with  a  study  of  the  moment  method  convergence  involving  the 
CAM-model  [25]. 

In  Quarter  3  for  CEM-model  it  was  obtained  the  analytical  solution  describing 
the  non-stationary  process  after  an  instant  switch-on  of  the  electric  field.  When 
obtaining  such  solutions,  the  interesting  perspectives  are  opened  for  the  convergence 
search.  Actually,  for  CEM-model  one  can  always  compare  the  DF  restored  via  the 
moments  with  the  exact  solution  and,  hence,  always,  an  opportunity  occurs  when 
the  moment  method  converges.  The  conclusions  obtained  for  CEM-model  turn  out 
to  be  true  for  other  interaction  potentials  of  ions  and  atoms  as  well,  a  problem 
concerning  the  DF  convergence  to  be  restored  via  the  moment,  analytically  built 
up,  being  separated  from  the  study  of  the  errors  which  arise  during  the  numerical 
solution  of  the  moment  system  itself. 

The  main  reason  of  the  moment  method  divergence  is  the  Grad  criterion  violation 
(7)  under  the  strong  deviation  from  equilibrium.  A  slow  convergence  takes  place 
under  rather  strong  deviations  from  equilibrium  when  the  system  is  the  nearest  to  a 
state  in  which  the  Grad  criterion  does  not  hold.  It  was  shown  analytically  [21]  that 
the  DF  at  the  stationary  state  (t  — >  oo)  does  not  satisfy  the  Grad  criterion  even  at 
a  low  value  of  the  electric  field  strength.  However,  in  a  course  of  the  non-stationary 
process  (at  finite  t),  the  DF  satisfies  the  Grad  criterion.  Clearly,  that  at  a  long 
time,  when  the  DF  differs  little  from  the  stationary  DF,  which  does  not  satisfy  the 
Grad  criterion,  the  number  of  the  moments  necessary  for  convergence  should  be  very 
increasing.  Any  question  concerning  the  true  state  of  art  can  be  answered  only  with 
the  relevant  calculation. 

When  solving  numerically  the  kinetic  problems,  the  cut-off  of  the  system  of  the 
moment  equations  is  necessary  to  be  carried  out,  i.e.,  instead  of  an  infinite  system 
of  the  moment  equations  the  finite  one  is  treated.  Cutting-off  the  moment  system, 
certain  domain  to  which  the  indices  l-  and  r  are  attributed  is  separated  out.  Under 
investigation  is  to  be  the  effect  of  the  form  of  such  a  domain.  The  aim  is  to  obtain 
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such  an  increment  of  the  moment  located  on  the  domain  boundary  being  dependent 
on  only  the  moments  located  in  its  interior.  In  this  case,  the  substitution  of  zero 
values  for  the  moments  in  its  exterior  has  not  any  influence  on  the  moment  increment 
within  this  domain. 

A  comparison  between  the  calculations  with  the  different  domain  of  summation 
was  carried  out.  Triangular  domain  REG  1 


l  <  N0  —  r 

(65) 

Triangular  domain  REG2 

r  <  (N0  -  1)12. 

(66) 

and  a  rectangular  domain  REG3 

r  <  Rq,  l  <  L0 

(67) 

Here,  Rq  is  the  number  of  the  Sonine  polynomials,  Lq  is  the  number  of  the  Legendre 
polynomials.  For  CEM-model,  the  most  natural  is  a  domain  REG2.  The  calculation 
results  are  not  inferior  than  those  with  a  choice  of  a  domain  REG  1  with  the  same 
amount  of  the  moments.  The  calculations  dealing  with  the  triangular  domain  are 
more  optimal  than  with  the  rectangular  one.  Nevertheless,  to  exclude  the  very  large 
values  of  /,  i.e.,  with  their  aid  to  decrease  further  the  number  of  the  moments,  the 
summation  can  be  carried  out  over  a  trapezoidal  domain  REGA 

l  <  N0  —  r,  l  <  Z/0,  L0  <  N0  (68) 

With  the  relevant  choice  of  N0  and  L0.  Further  calculations  involved  in  the  project, 
as  a  rule,  were  carried  out  with  a  summation  over  a  domain  REGA. 

Usually,  the  moment  method  is  applied  to  calculate  the  transport  coefficient. 

The  solution  of  the  system  of  the  moment  equations  with  the  comparatively 
large  number  of  the  moments  is  needed,  along  with  [4],  because  of  two  reasons: 
first,  the  larger  number  of  the  Sonine  polynomials  is  needed  to  be  considered  as 
compared  with  the  common  treatment  performed  via  the  simplest  variant  of  the 
Chapman- Enskog  method,  second,  when  increasing  the  electric  field  strength,  the 
senior  moments  can  affect  the  junior  ones  via  those  the  needed  transport  coefficients 
are  expressed. 

We  are  confronted  with  even  the  most  ambitious  problem,  i.e.,  to  build  up  the 
DF  with  the  moment  method.  Moreover,  we  propose  to  solve  this  problem  in  the 
non-  stationary  statement.  With  the  success  obtained,  the  transport  coefficients  will 
be  found  and  their  temporal  dependences  will  be  built  up.  Clearly,  the  calculation 
of  the  transport  coefficients  of  a  domain  in  the  plane  (/,  r)  only  will  be  of  essentially 
fewer  than  that  for  the  DF  calculation. 

To  build  up  the  DF  via  the  moments,  first,  fi(c)  (DF  in  /-subspace)  are  calcu¬ 
lated)  via  a  formula 

/I(c)  =  M(c)yC,Sr+1/2(c2)G,I  (69) 
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Then,  via  a  formula 

/(c)  =  J2Mc)pi(cos(e))  (70) 

the  DF  is  found. 

The  principal  studies  are  performed  at  the  axis  of  symmetry.  Consequently,  the 
numerical  studies  were  performed  for  the  cases  of  the  weak,  moderate  and  strong 
electric  field.  These  regions  are  determined  as  follows: 

Weak  electric  field  -  e  <  0.1, 

Moderate  electric  field  -  0.1  <  e  <  1, 

Strong  electric  field  -  e  >  1. 

Using  the  analytical  formulas  for  the  moments,  it  was  shown  that,  for  CEM- 
model,  the  moment  method  results  in  a  divergent  series  in  attempting  to  build  up 
the  DF  at  any  /  and  any  value  of  the  electric  field  e. 

However,  at  any  finite  time,  the  absolute  values  of  the  moments,  i.e.,|CV/)  de¬ 
crease  with  an  increase  in  r  at  any  /,  and  the  DF  expansion  in  terms  of  the  moments 
should  principally  be  convergent. 

When  investigating  the  weak  electric  field,  the  problem  of  the  DF  building-up  in 
the  stationary  state  receives  primary  emphasize.  It  turns  out  that,  in  this  region  of 
e,  there  is  the  minimum  on  the  | Crj\(r)  dependences.  Even  on  the  upper  boundary 
of  this  region  of  e  =  0.1,  the  minimum  is  very  deep,  at  the  minimum,  [C^oD  is  12 
orders  lesser  than  | C'o.o  | ;  the  minimum  coordinate  r  =  25. 

At  the  small  values  of  e,  a  departure  from  equilibrium  is  small.  Hence,  the 
number  of  the  moments  needed  to  build  up  the  DF  should  be  small.  It  results  in 
a  real  possibility  to  build  up  the  stationary  DF  via  the  moments  involving  the  true 
choice  of  the  cut-off  domain.  As  one  would  expect,  the  cut-off  should  be  performed 
somewhere  within  the  minimum  neighbourhood. 

When  transferring  in  a  domain  of  the  moderate  field,  with  an  increase  in  e, 
there  is  a  very  fast  decrease  in  a  minimum  depth  on  a  \Cr)i\(r)  dependence.  An 
attempt  to  build  up  the  stationary  DT  via  the  stationary  moment  does  not  result  in 
a  satisfactory  yield  already  at  e  =  0.25.  In  the  iteration  method,  advanced  by  Kihara 
[22],  to  built  up  the  stationary  solution,  the  moments  need  to  be  found.  Emphasize, 
that,  for  CEM-model,  the  stationary  moments  are  known.  Nevertheless,  there  is 
not  a  grain  of  hope  to  succeed  in  building-up  the  stationary  DF  in  a  domain  of  the 
moderate  field  via  the  iteration  method.  This  does  not  ruled  out  the  possibility  to 
find  several  moments  with  the  small  indices,  particularly,  a  drift  velocity  with  which 
the  mobility  can  be  found  using  this  method. 

Quite  to  the  contrary  is  a  deal  concerning  the  DF  restoration  via  the  moments 
in  a  course  of  non- stationary  process.  It  is  due  to  a  form  of  \Crj\(r)  dependences. 
In  any  finite  time,  beginning  with  certain  r,  the  moments  decrease.  First,  there  is 
observed  the  minimum  at  |CV/(r)  dependence,  then,  there  is  the  maximum.  When 
increasing  t,  the  coordinate  of  the  maximum  and  its  height  increase  with  no  limit. 
At  the  maximum  left,  there  is  a  region  r,  in  which  interior  a  dependence  Crj(r)% 
corresponding  to  the  non-stationary  solution,  coincides  with  the  stationary  one,  and 
this  region  is  enlarged  with  an  increase  in  t. 
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For  the  convergence  of  a  series  (69)  to  the  exact  solution,  one  needs  to  pass  the 
maximum,  when  summing  over  r,  and  enter  a  region  r  where  the  values  of  \Crj \ 
are  sufficiently  small.  At  very  large  time  when  the  coordinates  of  the  maximum  are 
large,  such  a  summation  turns  out  to  be  very  complicated.  However,  in  the  domain 
of  the  moderate  electric  held,  the  DF  restoration  is  successfully  performed  up  to 
sufficiently  long  time. 

In  [21],  it  was  shown  that  the  process  to  maintain  the  stationary  state  is  ac¬ 
companied  by  propagation  of  a  front  along  c^-axis  in  the  velocity  space.  Behind 
this  front,  the  stationary  state  develops  almost  immediately.  In  a  domain  of  the 
negative  velocities,  the  stationary  state  yielding  is  over  approximately  by  t  —  1.  By 
the  sufficiently  long  time,  the  DF  coincides  with  the  stationary  one  within  a  very 
wide  region  cz. 

As  a  result,  with  the  aid  of  ”non-stationary  moment  method”,  we  succeed  in 
solving  the  seemingly  unsolvable  problem  to  build  up  the  stationary  solution  in  a 
wide  region  of  cz. 

In  a  domain  of  the  strong  electric  held,  the  problem  of  passing  the  maximum 
turns  out  to  be  more  complicated.  For  the  convergence,  the  cut-off  region  should  be 
enlarged  strongly  in  a  (r,  /)-plane. 

When  calculating  fi(c),  there  is  a  product  of  the  Sonine  polynomials  Sl+l,2(c2) 
and  the  moments  Crj  under  a  sign  of  a  sum  (69).  The  moments  have  the  alternating 
sign  with  a  change  in  r,  this  is  true  for  the  Sonine  polynomials  as  well.  As  a  result, 
the  terms  of  the  sum  change  their  signs  and  the  sum  turn  out  to  be  many  orders 
lesser  that  the  absolute  value  of  a  maximum  term  of  the  sum.  If,  when  passing  the 
maximum,  a  ratio  of  the  maximum  general  term  to  the  exact  value  is  greater  than 
1016,  the  loss  in  accuracy  occurs  and  the  true  value  is  out  of  reach,  the  number  of 
the  terms  of  the  sum  being  made  as  large  as  is  wished. 

It  was  maintained  that,  if  e  and  t  are  simultaneously  changed  but  in  the  manner 
of  a  product  e  •  t  being  constant,  the  dependences  \Crj\(r )  turn  out  to  be  very 
similar.  Thus,  certain  parameter  of  similarity  is  found  -  e  •  t.  It  is  shown  that  a 
good  convergence  takes  place  at  e  •  t  <  5. 

In  a  domain  of  the  strong  electric  held,  the  ”  non-stationary  moment  method”  to 
build  up  the  stationary  regime  is  improved. 

The  analytical  solutions  for  CEM-model,  being  built  up,  simultaneously,  both 
for  the  Boltzmann  equation  itself  and  the  moment  system,  gave  an  insight  into 
capabilities  of  the  standard  moment  method. 

In  Fig.  2,  for  CEM-model,  the  DF  at  the  axis  of  symmetry  (cp  =  0)  is  shown 
for  different  times  at  e  =  2.  It  shows  in  what  manner  increases  the  number  of  the 
Sonine  Rq  and  Legendre  L0  polynomials  to  be  taken  into  consideration  in  the  DF 
expansion. 

For  cz  >  —4,  there  is  a  fine  coincidence  with  the  corresponding  analytical  solu¬ 
tions.  For  cz  <  —4,  a  very  fast  yield  toward  the  stationary  solution  but,  as  a  time 
increases,  the  more  difficulties  arise  when  building  up  the  solution  via  the  moment 
method. 
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Figure  2:  The  distribution  function  at  the  symmetry  axis  at  different  times  for 
CEM-model  at  e  =  2.0. a-  f(cz,t).  b  -  f(cz,t)/f(cz,  oo).  Curves  1-5  refer  to  a  time 
t  =  0.5  (I/0  =  35,  Rq  =  23);  1.5  (L0  =  30,  R0  =  26);  2.5  (L0  =  37,  R0  =  70);  3 
(L0  =  75,  Rq  =  100);  5  (L0  =  100,  i?0  =  280). 

Thus,  with  an  increase  in  the  electric  field,  the  DF  building-up  via  the  moment 
method  is  limited  by  the  shorter  times,  and,  to  solve  the  problems  at  the  very  strong 
fields,  it  is  necessary  to  develop  the  new  methods  in  kinetic  problems  solution.  In 
spite  of  the  above  objection,  basing  on  the  studies  performed,  the  conclusions  can 
be  derived  as  follows:  1.  The  non-stationary  moment  method  to  build  up  the  DF 
is  developed.  This  method  is  an  advance  in  the  standard  moment  method  involving 
the  DF  expansion  about  the  Maxwellian  of  which  temperature  coincides  with  the 
temperature  of  the  background  atoms.  2.  Using  the  non-stationary  moment  method 
,  we  succeed  in  building  up  the  solution  in  a  domain  of  the  moderate  and  strong 
electric  fields.  3.  At  least,  by  one  order  of  value  is  increased  e  for  which,  using  the 
moment  method,  the  DF  turns  out  to  be  build  up  at  the  stationary  state. 

5  Fifth  quarter 

In  Quarter  5  [26],  the  solutions  of  the  non-stationary  problem  were  obtained  under 
a  sharp  switch-on  of  the  electric  field  for  the  ions  with  several  interaction  models. 
In  the  detailed  5th  quarterly  report,  there  are  many  detailed  proofs  and  research 
results  as  well  as  many  illustrations.  The  Quarter  5  being  terminated,  this  detailed 
report  was  sent  to  EOARD  as  well  as  the  brief  report. 

In  the  first  section  of  the  5th  quarterly  report  [26],  with  an  example  of  CEM- 
model,  the  studies  on  the  moment  method  convergence  were  continued. 

In  [25],  it  was  revealed  that,  in  the  case  CEM-model,  an  increase  in  the  elec¬ 
tric  field  e  as  well  as  a  time  t  result  in  poorer  convergence  of  the  DF  restoration 
via  the  moments.  It  is  because  of  the  high  maximum  arising  on  a  dependence  of 
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the  moments  \Cri\  on  r  at  a  fixed  /;  it’s  height  and  it’s  r-coordinate  increase  in 
time.  When  calculating  the  DF,  one  should  go  essentially  to  the  right  of  this  max¬ 
imum,  The  method  of  the  MEs  calculation  of  large  indices  (matrix  elements  of  the 
collision  integral)  developed  by  us  is  ready  to  make  it.  However,  if  the  degree  of 
non-equilibrium  is  high  and  the  moment  values  of  the  maximum  are  enormous,  it 
results  in  a  loss  in  accuracy  of  the  DF  calculation  during  summation  as  well  as  in  the 
convergence  default.  It  occurs  when  the  DF  is  slightly  distinct  from  the  stationary 
one  which  does  not  satisfy  the  Grad  criterion.  In  the  5th  quarterly  report,  the  study 
of  the  attempts  to  ameliorate  the  moment  method  convergence  clue  to  a  change  in 
basis,  i.e.,  an  increase  in  the  temperature  of  the  weight  Maxwellian,  were  studied. 
Earlier,  the  need  of  the  basis  change  was  mentioned  in  [5],  [9],  [10].  However,  in 
these  works,  the  ME  building-up  in  the  new  basis  was  performed  via  the  calculation 
of  the  complicated  integrals  with  the  large  time  consumption  [9]  with  no  certainty  in 
a  good  accuracy,  the  ME  building-up  being  successfully  carried  out  with  the  rather 
small  indices  (/  <  7,  r  <  10)  only. 

When  investigating  the  matrix  elements  of  the  collision  integral  in  [1] ,  a  point  on 
a  transfer  from  one  basis  to  another  was  considered  in  details.  The  transfer  matrix 
(operator  D)  was  built  up,  the  bases  being  distinct,  in  the  general  case,  not  only 
in  the  temperature  but  the  mean  velocity  of  the  Maxwellian.  When  building  up  an 
operator  D ,  we  use  a- u-r epr esent at io n  of  the  Boltzmann  equation  [17],  [27].  In  the 
report,  the  detailed  derivation  of  a  matrix  D,  which  relates  both  bases,  which  are 
distinct  only  in  their  temperature.  With  the  matrix  D,  one  can  obtain  the  ME  in 
the  new  basis  via  the  ME  in  the  initial  basis.  Using  the  results  of  the  2nd  quarter, 
we  can  pass  from  the  global  basis  to  the  local  one  and  solve  the  moment  system  in 
the  new  basis;  we  call  such  approach  the  modified  moment  method.  In  general, 
the  parameters  of  the  weight  Maxwellian  experience  the  changes  all  the  time  in  a 
course  of  the  non-  stationary  process.  To  verify  the  capabilities  of  this  method,  there 
is  no  need  in  solving  the  system  of  the  moment  equations  in  the  new  basis  with  the 
temporal  temperature  dependence.  One  can  transform  the  moments,  which  were 
built  up  in  an  initial  basis  C°h  to  the  new  basis  and  obtain  . 

Thus,  when  transferring  from  a  basis  T0  to  the  T)  one,  the  moments  change  as 
follows 

OO 

Cl,=  Y.DrlM(TxA)Ci„  (71) 

k= 0 

It  was  considered  an  example:  e  =  1,  t  —  5.  With  such  a  field  and  time,  the 
moments  Cri  in  a  basis  T0  are  very  large. 

It  was  developed  a  special  code  in  which  the  moments  built  up  in  a  basis  Tq 
(BasO)  were  transformed  via  a  formula  (71)  to  a  basis  Tf  (Basl).  The  relevant 
parameter  for  this  transformation  is  a  ratio  RT  =  T\  /T0. 

In  Fig.  3,  the  calculation  results  at  RT  =  1.5  are  presented. 

One  can  see  in  this  figure  that,  in  a  subspace  l  =  0,  a  transfer  to  the  new  basis  is 
accompanied  with  a  change  in  the  Cri  maximum  by  ~  5  orders  of  value,  in  a  subspace 
l  =  10,  by  ~  3  orders.  In  both  subspaces,  there  is  a  shift  of  the  distribution  Cri 
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toward  the  domain  of  small  r  —  a  coordinate  of  the  maximum  is  strongly  decreases, 
i.e. ,  at  l  =  0,  it  changes  from  24  to  9.  Thus,  even  a  comparatively  small  change  in 
the  basis  temperature  results  in  an  essential  facilitation  in  the  DF  restoration  via 
its  moments. 


Figure  3:  The  dependence  of  the  moments  Cri  on  r  for  CEM-model  in  two  bases 
with  the  temperatures  T0  and  Ti;  RT  =  1.5,  L0  =  37,  Rq  =  70.  Figure  a  -  BasO, 
l  —  0,  b  -  Basl,  l  —  0,  c  -  BasO,  l  =  10,  d  -  Basl,  l  =  10. 


To  build  up  the  DF  with  the  same  accuracy  as  for  the  initial  basis,  it  is  needed 
essentially  smaller  number  of  the  moments.  The  results  obtained  lend  support  to 
the  utility  in  involving  the  modified  moment  method. 

Our  experience  in  the  moment  system  transformation  from  one  basis  to  another 
gives  the  way  to  solve  this  problem  by  the  efforts  of  our  collective  in  the  case  of  a 
new  project  devoted,  generally,  to  the  development  of  the  modified  moment  method. 


In  the  5th  quarter,  also,  the  formulas  relating  the  physical  moments  and  the 
moments  Cri  were  obtained.  Beside  of  such  moments  as  a  drift  velocity  Ud  and 
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mobility  K  =  ud/e,  the  expressions  for  a  total  energy,  total  energy  flow,  thermal 
energy,  and  heat  flow. 

In  the  axially  symmetric  case,  an  expansion  of  the  angular  dependence  of  the 
DF  is  carried  out  only  in  terms  of  the  Legendre  polynomials. 

The  first  physical  moment,  i.e.,  concentration  coincides  with  Coo-  In  the  process 
under  consideration,  n  experiences  no  changes,  and  we  set  C0 o  =  1,  i.e.,  the  DF  is 
said  to  be  normalized  by  unity.  The  physical  moment  are  the  mean  values  of  the 
physical  quantities  which  are  expressed,  usually,  via  the  different  combinations  of 
the  various  powers  of  the  particle  velocities.  Mean  values  to  be  normalized  by  unity 
are  determined  as  follows 

(444)  =  /  /444<*3e-  (72) 

A  mean  velocity  of  the  ions  is  their  drift  velocity  to  be  denoted  by  ud.  For  it  we 
have 

ud  =  (c2)  =  ^0)1  (73) 

Determine  the  mean  energy  of  the  ions  W  =  (mv2/ 2)  in  a  frame  related  to  the 
atoms. 

For  dimesionlcss  energy  we  have 

^  =  3/2 (Coo  -  C10)  (74) 

The  mean  dimensionless  energy  flow  in  the  frame  is  expressed  via  Cri  as  follows 

- t2_  =  5/4(670!  -  Cn)  (75) 

kT^2kT  /m 

To  determine  the  mean  thermal  energy,  mean  heat  flow  and  other  moments 
related  with  a  chaotic  motion  of  the  particles,  it  is  convenient  to  involve  a  ”  proper” 
velocity  of  the  particles 

C  =  c  -  (c)  =  c  —  ud  (76) 

The  ions  thermal  energy  —  the  energy  in  a  frame  related  to  the  ions  is  determined 
as  follows 

WT 

W  = {c)'  (77) 

and  a  heat  flow 

Qt/ (kT \J  2kT /m)  =  ( C2CZ )  (78) 

After  a  series  of  transformations  we  obtain 

|^=3/2(Cm-C1„)-l/4C021,  (79) 

Qr/{kT  J2kr]m)  =  -5/4C„  +  5/4 C10C01  +  1/4C03,  -  1/2 C„2C01  (80) 
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T77T  —  0)0  —  Cio  —  I/2C02  (84) 

kl 

The  physical  moments  behaviour  is  calculated  after  a  sharp  switch-on  of  the 
electric  field  for  several  interaction  models: 

Pseudo-Maxwellian  molecules,  a  total  cross  section  is  inversely  proportional  to 
the  relative  velocity  E  =  C/g,  angular  part  of  the  cross  section  is  isotropic 

CEM-model,  the  total  cross  section  is  also  inversely  proportional  to  the  relative 
velocity,  angular  part  of  the  cross  section  —  the  resonance  charge  exchange  refers 
to  180°-  scattering. 

CEHS-model,  the  total  cross  section  does  not  depend  on  a  velocity  E  =  const , 
angular  part  of  the  cross  section  is  the  same  as  for  CEM-model  of  the  charge  ex¬ 
change. 

HS-model,  hard  sphere  model,  the  total  cross  section  does  not  depend  on  a 
velocity  E  =  const,  angular  part  of  the  cross  section  -  isotropic  scattering. 

Such  normalization  is  advanced  that,  in  the  limit  of  very  weak  fields,  the  mobility 
equals  unity.  Such  normalization  is  named,  further,  non-standard  normalization. 


ft  was  shown  that,  with  the  diagonal  interaction  matrix  (CEM-model  and  pseudo- 
Maxwellian  molecules),  the  temporal  dependence  of  the  ion  mobility  after  the  field 
switch-on  is  the  universal  curve  and  does  not  depend  on  e  as  well  as  the  limit  (sta¬ 
tionary)  value  of  mobility.  For  other  two  models  this  dependence  changes  with  an 
increase  in  e.  The  limit  (stationary)  value  of  mobility  depends,  also,  on  the  electric 
field,  i.e. ,  it  decreases  with  an  increase  in  e. 

The  transformation  to  the  non-standard  normalization  occurs  as  follows:  the 
problem  is  solved  using  the  standard  normalization  at  very  small  e,  the  mobility 
is  found,  being,  in  general,  non  unity,  this  value  is  a  coefficient  by  which  all  the 
MEs  are  divided,  Clearly,  that,  changing  the  MEs,  we  reform  the  coefficient  in  the 
collision  integral,  i.e.,  we  change  the  time  unit  r.  As  a  result,  both  the  time  scale 
and  the  value  of  the  dimensionless  electric  field  e  change.  After  the  transfer  to 
the  non-standard  normalization,  for  all  interaction  models,  the  ion  mobility  goes  to 
unity,  K  — >  1  at  e  — >  0. 
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Figure  4:  The  dependence  of  the  standard  mobility  K  on  a  strength  e  for  HS-model. 
For  comparison,  the  dashed  curve  shows  analogius  dependences  for  CEHS-modcl. 
The  non-standard  normalization  is  used. 

If  for  the  models  with  the  constant  collision  frequency  the  mobility  does  not 
depend  on  e  and  equal  always  unity  in  the  non-standard  normalization,  for  the 
models  with  A  =  const,  the  mobility  decreases  with  an  increase  in  e,  such  decrease 
occurring  more  rapidly  than  for  the  hard-sphere  model,  as  one  can  see  in  Fig.  4,  in 
the  case  of  CEHS-model. 

In  a  course  of  investigations,  we  found  a  series  of  the  interesting  details  concerning 
the  behaviour  of  a  heat  flow  and  transverse  component  of  thermal  energy,  and  it  was 
shown  that  their  behaviour  changes  essentially  when  transferring  from  one  model  of 
interaction  to  the  other  one. 

Thus, 

1.  The  comprehensive  investigations  of  the  moment  method  convergence  are 
finished  for  CEM-model. 

2.  The  way  to  ameliorate  the  moment  method  convergence  are  outlined  involving 
the  transfer  to  the  modified  moment  method  in  the  cases  when  the  DF  does  not 
satisfy  the  Grad  criterion. 

3.  For  several  interaction  models,  the  DF  calculations  were  carried  out  as  well  as 
its  moments  in  a  domain  of  the  moderate  electric  fields  as  well  as  at  the  beginning 
of  a  region  of  the  strong  ones. 


6  Sixth  quarter 

During  the  6th  Quarter  [28]  our  work  was  carried  out  concerning  two  items. 

First,  the  study  of  the  ion  DF  under  the  electric  held  was  continued  for  the 
various  models  of  interaction. 
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Second,  a  system  of  the  moment  equations  was  deduced  for  the  ions  moving  in 
the  crossed  electric  and  magnetic  fields.  In  this  case,  the  DF  is  no  longer  the  axially 
symmetric  one,  and  the  equations  for  the  most  general  case,  when  the  DF  is  the 
3-Dimension  one  in  a  velocity  space,  should  be  used. 

Namely,  in  the  first  part  of  the  6th  quarterly  report,  the  ultimate  choice  of  the 
standard  normalization  was  made  about  which  we  wrote  above  when  describing  the 
Quarter  3.  For  the  models  with  A  =  const,  the  definition  of  the  dimensionless 
electric  held  was  made  in  a  form  (54). 

Under  investigation  is  the  DF  evolution  at  various  values  of  e  and  different  inter¬ 
action  models.  The  domains  of  the  weak,  moderate,  and  strong  electric  held,  where 
the  Chapman-Enskog  method  [29]  is  applicable,  were  considered  in  consecutive  or¬ 
der.  In  these  cases,  there  are  many  known  results  ready  for  comparison. 

At  the  very  small  e,  in  the  Chapman-Enskog  method,  the  DF  is  sought  in  a  form 

/(c)  =  Af(c)(l  +  e<p(c)).  (85) 

Here,  M(c)  is  the  Maxwellian.  We  build  up  the  DF  via  moment  method  involving 
very  large  number  of  the  MEs  of  the  collision  integral.  It  gives  an  opportunity  to 
build  up  the  DF  with  a  high  accuracy,  particularly,  for  the  very  small  values  of  e. 
Always,  so  small  e  can  be  found  that  a  correction  ip  does  not  depend  on  e.  This  tp, 
built  up  in  such  a  manner,  we  call  the  universal  Chapman-Enskog  correction. 

We  calculated  the  temporal  dependences  of  the  universal  ip  for  four  models.  If,  for 
the  CEM-model  and  the  pseudo-Maxwellian  molecules  ,  this  function  is  proportional 
to  c  (the  coefficient  of  proportionality  depends  on  a  time),  a  tendency  to  saturation 
is  inherent  for  CEHS-  and  HS-models  at  longer  times  and  high  values  of  cz.  The 
value  of  this  saturation  is  about  1.5  times  lesser  for  CEHS-model  as  compared  with 
HS-model. 

When  studying  the  processes  under  the  very  weak  held,  it  was  shown  that  the 
DF  building-up  requires  the  DF  expansion  with  the  large  number  of  the  Sonine 
polynomials  (N0  =  128  and  more)  although  the  ion  mobility  is  easily  determined 
with  the  small  amount  of  the  moments  (2  or  3). 

There  is  the  large  section  in  the  6th  quarterly  report  devoted  to  another  method 
in  which  there  are  no  principal  limits  of  the  moment  method  related  with  the  Grad 
criterion.  Note,  that  the  Grad  criterion  arises  at  the  stage  of  an  expansion  in  terms 
of  the  Sonine  polynomials  which  are  orthogonal  with  the  Maxwellian  weight.  We 
suppose  to  carry  out  the  expansion  in  terms  of  the  spherical  harmonics  (in  the  case 
under  consideration,  in  terms  of  the  Legendre  polynomials).  This  method  of  solution 
of  the  Boltzmann  equation  goes  back  to  the  works  [30],  [31]. 

In  the  axially  symmetric  case,  the  DF  is  represented  as  (70).  At  very  small 
e,  when  the  representation  (85)  is  true,  a  function  (p  is  expanded  in  terms  of  the 
Legendre  polynomials,  only  <p\  being  non-zero.  In  [30],  [31],  the  linearized  Boltz¬ 
mann  equation  was  considered  for  HS-model.  It  was  shown  that  outgoing  term  of 
the  collision  integral  corresponds  to  multiplication  of  the  DF  by  a  function  k(c )  at 
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a  scattering  on  equilibrium  background  gas.  Hilbert  found  a  function  k(c),  in  the 
standard  normalization  it  is  as  follows 

_c2 

Kc)  =  ^172  +  (l/2c  +  c)$(c).  (86) 

Incoming  term  of  the  collision  integral  takes  a  form  of  an  integral  of  the  DF  with 
a  linear  kernel,  of  which  explicit  expression  was  obtained  in  [30]  and  the  expansion 
in  terms  of  the  spherical  harmonics  in  [31].  Using  these  results  for  small  e,  the 
Boltzmann  equation  can  be  easily  written  as  follows 

-2 c  =  J kf(c,  ci)</?i(ci,  t)c\dci  -  k{c)ipx{c,t).  (87) 

0 

For  CEM-model,  kf(c,ci)  =  0,  and,  k(c)  =  lat  the  stationary  state,  we  have 

V?i(c)  =  2c,  (88) 

that  coincides  with  the  numerical  solution. 

If  a  perturb  of  the  DF  from  the  equilibrium  one  is  not  small  (in  the  case  under 
consideration,  e  is  not  small),  the  Boltzmann  equation  is  substituted  by  the  system 
of  equations  for  fi(c).  For  the  linear  Boltzmann  equation  (a  small  admixture  on  the 
equilibrium  background),  we  have 

1  l^1  =  J {c,  ci)fi(ci,t)cldci  -  k(c)fi(c,t).  (89) 

0 

Here,  Di/Dt  corresponds  to  the  left  side  of  the  Boltzmann  equation,  and  the  kernels 
in  (87)  and  (89)  are  interrelated  via  a  simple  relationship 

K  (c,  Cl)  =  - (c,  d)M (ci).  (90) 

A  great  advantage  of  the  transfer  from  the  system  of  the  moment  equations  to 
the  system  (89),  besides  the  removal  of  limitations  related  with  the  grad  criterion, 
is  the  transfer  in  the  outgoing  term  to  a  very  simple  multiplication  operator  by  a 
function  k(c).  In  the  moment  method,  the  outgoing  term  contribution  prevails. 

In  the  case  of  HS-model,  the  kernels  kf(c,  ci)  and  kf(c,  ci)  are  easily  expressed 
via  the  symmetrical  kernels  obtained  by  Hecke  [31].  However,  for  the  arbitrary  cross 
sections  of  interactions,  there  are  no  analytical  expression  for  the  kernels. 

In  [1],  it  was  shown  that  the  kernels  can  be  build  up  via  the  known  MEs.  The 
recurrence  relationships  with  which  aid  we  build  up  the  ME  [3]  are  true  for  the  ME 
of  the  incoming  term  of  the  collision  integral.  In  [1],  it  was  shown  that  the  kernels 
can  be  built  up  for  arbitrary  cross  sections,  even  for  the  non-  linear  Boltzmann 
equation. 


Figure  5:  The  DF  at  different  times  for  the  hard-sphere  model  at  e  =  0.5  in  the 
non-standard  normalization.  L0  =  64,  N0  =  128.  Figure  a  -  f(cz,t),  Fig.  b  - 
f(cz,  t)/f(cz,  oo).  The  curves:  1  -  t  —  0.1,  2  -  t  =  0.5,  3  t  =  1,  4  t  =  3,  5 
t  —  5,  6  -  t  —  10.  Dashed  curve  refers  to  t  —  0. 


In  the  linear  case  of  our  interest,  we  have 


OO 

fcz+(c,ci)  =  M(d)clc[  J2 

r, ri=0 


,1 


(91) 


where  arj  is  a  square  of  a  norm  of  the  Sonine  polynomial. 

Thus,  when  developing  this  way,  there  are  actual  possibilities  of  an  insight  into 
domain  of  the  very  strong  electric  field. 

Using  a  kernel  for  the  hard-sphere  model,  the  estimations  were  carried  out  and 
the  asymptotic  value  of  <^i(c)  was  found  at  c  — >  oo.  In  the  standard  normalization, 
it  turns  out  to  be  equal  to  6. 

Further,  in  the  first  part  of  the  6th  quarterly  report,  using  the  standard  moment 
method  for  several  models  of  interaction,  the  DF  was  under  study  in  a  domain  of 
the  moderate  and  strong  electric  field. 

As  the  examples,  there  are  two  figures.  In  Fig.  5a j  is  presented  the  DF  evolution 
for  HS-model  at  e  =  0.5,  in  Fig.  5b,  the  same  evolution  is  shown  involving  the 
relative  DF  where  f(cz)  refers  to  the  stationary  DF,  and  as  the  stationary,  f(cz)  is 
chosen  at  t  —  10. 

In  Fig.  6,  the  stationary  DFs  are  shown  for  a  series  of  e  (HS-model). 

The  main  results  of  this  part  can  be  formulated  as  follows 

1.  The  DF  was  built  up  successfully  via  the  standard  method  in  a  wide  region 
of  e. 

2.  An  idea  to  solve  the  non-standard  problem  at  the  instant  switch-on  of  the 
electric  field  turns  out  to  be  the  successful  one.  Even  at  the  strong  field,  when 
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Figure  6:  The  DF  f(cz,  oo)  at  different  e  for  the  hard-sphere  model  in  the  non¬ 
standard  normalization.  Lq  =  128,  Nq  =  128.  The  curves:  1  -  e  =  0.2,  2  -  e  =  0.5, 
3  -  e  =  1,  4  -  e  —  1.5,  5  -  e  =  2.  Dashed  curve  refers  to  t  —  0. 

the  stationary  solution  does  not  satisfy  the  Grad  criterion  and  can  not  be  built  up 
via  the  solution  of  the  standard  moment  system,  we  succeed  to  nearly  reach  the 
stationary  solution. 

3.  For  the  first  time,  it  was  shown  in  what  manner  the  DF  evolution  proceeds 
in  a  course  of  the  non-stationary  process  and  how  it  is  affected  by  the  interaction 
model. 

4.  The  standard  moment  method  is  far  from  being  exhausted,  and  can  be  used 
for  solving  the  more  complicated  problems,  namely,  concerning  the  ion  kinetics  in 
crossed  electric  and  magnetic  fields. 

5.  The  ways  to  overcome  those  limits  which  are  relevant  for  the  standard  mo¬ 
ment  method  ere  outlined.  In  the  second  part  of  the  6th  quarterly  report,  under 

consideration  are  the  moment  equations  describing  the  ion  kinetics  under  crossed 
electric  and  magnetic  fields  when  the  DF  expansion  should  be  performed  over  a  full 
set  of  the  spherical  Hermite  polynomials  (2). 

From  the  collision  integral  invariance  relative  to  rotation,  the  relationships  be¬ 
tween  the  MEs  of  different  indices  m  and  i  were  derived.  It  was  shown  that  the 
MEs  with  arbitrary  indices  m  and  i  are  proportional  to  the  corresponding  axially 
symmetric  MEs  (m  =  0,  i  =  0)  with  the  coefficients  of  proportionality  expressed  via 
the  Klcbsch-Gordan  coefficients. 

The  main  result  for  the  case  under  our  consideration  can  be  formulated  as  follows. 
The  linear  matrix  element  -f^qmiooo  is  non  zero  only  at  li  =  /,  mi  =  m, 
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and  this  ME  equals  the  corresponding  axially  symmetric  ME: 


jy-rlm  _  7^rZ 

■^■ri/mOOO  -^riZOO 


=  A 


r,ri,Z* 


(92) 


Thus,  the  number  of  the  axially  symmetric  MEs  calculated  earlier  is  well  suit  for 
solving  the  problem  concerning  the  ion  motion  under  crossed  electric  and  magnetic 
fields.  Only  the  number  of  the  equations  increases,  i.e.,  at  each  /,  an  index  m  changes 
from  0  to  l. 

Under  crossed  electric  and  magnetic  fields  the  Boltzmann  equation  in  the  dimen¬ 
sion  form  is  as  follows 


d£ 

dt 


v-|  +  -(E  +  vxH) 


dr 


m 


d£ 

dv 


=  J, 


col • 


(93) 


Under  consideration  are  the  spatially  uniform  case  and  the  mutually  perpendicular 
electric  and  magnetic  fields;  for  concreteness,  the  electric  field  is  directed  along  ,2- 
axis  and  the  magnetic  one  along  y- axis.  Using  the  standard  notion  of  the  unit 
vectors  along  these  axes:  k  and  j,  we  have  in  the  dimensionless  form 


_L  +  (ek  +  /,(cxj))._L  =  jcoh  (94) 

Here,  the  dimensionless  magnetic  field  h  is  proportional  to  H 

6  6 

h  =  — Hr,  h  =  lult,  <^l  =  — H ,  (95) 

m  m 

where  ujl  is  a  Larmor  frequency.  Remind  that  an  expansion  of  the  angular  depen¬ 
dence  of  the  DF  is  carried  out  in  terms  of  the  real  spherical  harmonics  Y}lm : 

Yi°m(0,  <p)  =  Pim(x) cos  rmp,  Y^n(d,  <p)  =  P/n(x)sin  imp,  x  =  cos  9  ,  (96) 


and  Pim(x)  are  the  associated  Legendre  polynomials.  Here,  an  angle  6  is  counted 
from  z-axis,  and  an  angle  <p  -  from  a  ( x ,  2:) -plane. 

In  the  2nd  quarterly  report  [19] ,  the  left  sides  of  the  moment  equations  were  built 
up  when  supposing  the  dependence  of  the  external  field  on  a  velocity.  To  obtain  an 
expression  for  the  left  side  of  the  moment  equations  in  the  case  under  consideration, 
the  rather  cumbersome  calculation  needs  to  be  carried  out. 

As  a  result,  the  system  of  the  moment  equations  was  built  up  describing  a  be¬ 
haviour  of  the  charged  particles  under  crossed  electric  and  magnetic  fields,  ft  takes 
a  form 

h 

-g*  +  2  ((*  -  m W  +  m+  l)Crlm+ 1  -  dm-lC^P 
e^2l  +  3  r(^  +  +  _  2/  —  1^  _  ^ r,ri,lClriim ,  (97) 

dm  1  ^mO  • 
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7  Seventh  quarter 

The  main  problem  to  be  solved  in  the  7th  quarter  is  to  build  up  the  analytical 
solutions  for  several  first  moments  in  crossed  electric  and  magnetic  fields  [32], 

Under  investigation  is  a  transient  process  after  the  switch-on  of  the  electric  field 
under  existing  magnetic  field  being  perpendicular  to  the  electric  one. 

At  the  symmetric  initial  conditions  under  our  consideration,  an  index  i  in  a 
system  (97)  can  be  only  zero  and,  usually,  it  is  omitted. 

Beside  of  the  analytical  solutions,  a  code  for  numerical  solution  was  developed 
(97).  The  analytical  solutions  for  the  moments  Crim  were  successfully  built  up  for 
the  molecules  of  the  Maxwellian  type,  i.e.,  those  of  which  cross  section  is  inversely 
proportional  to  the  relative  velocity  and  its  angular  dependence  is  the  arbitrary 
one.  The  non-standard  normalization  is  used,  all  the  MEs  being  divided  by  the 
normalization  coefficient  which  is  chosen  to  satisfy  the  condition  when  the  mobility 
equals  unity  and  coincides  with  that  for  CEM-model  at  zero  magnetic  field  and  very 
weak  electric  field.  For  the  Maxwellian  molecules  with  arbitrary  angular  part  of  the 
scattering  cross  section,  the  matrix  corresponding  to  the  linear  collision  integral  is 
diagonal,  all  diagonal  MEs  are  negative  and  do  not  depend  on  an  index  m.  For  such 
molecules,  the  renormalization  coefficient  turns  out  to  be  equal  to  an  initial  Aooi, 
i.e.,  in  the  new  normalization,  Aq01  =  —1.  We  introduce  a  new  notion 

A  ri  =  -A  *  .  (98) 

For  the  arbitrary  Maxwellian  molecules  in  the  non-standard  normalization,  Aoi  =  1. 

Note,  that,  with  a  change  in  the  normalization,  the  time  and  electric  field  units 
change  also  change  as  well  as  that  of  the  magnetic  field  in  the  case  under  consider¬ 
ation. 

For  concreteness,  the  pseudo-Maxwellian  molecules  were  considered  with  the 
isotropic  angular  part  of  the  scattering  cross  section. 

If,  with  no  magnetic  field,  a  drift  velocity 

u d  =  f  fed3  c 

has  a  single  component  Udz  related  with  the  moment  Coio  via  a  formula  Udz  =  u\\  = 
Coio/2,  at  the  mutually  perpendicular  electric  and  magnetic  fields,  another  compo¬ 
nent  of  a  drift  velocity  vector  Udx  arises  directed  along  x-axis  being  perpendicular  to 
the  directions  of  the  electric  and  magnetic  fields  and  equals  to  Udx  =  u±  —  C0n/ 2. 
Namely,  a  component  of  the  drift  velocity  Udx  determines  a  Hall  current.  Besides  the 
drift  velocity  vector,  a  mobility  vector  K:  Udz  =  AT,  Udx  =  Kh.aie  can  be  introduced. 

Usually,  to  characterize  a  direction  of  a  drift  velocity,  a  so  called  Lorentz  angle 
( p  is  introduced  which  is  determined  as  follows 

tan  ip  =  — .  (99) 

U\\ 
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The  value  of  the  drift  velocity  is  determined  as 


(100) 


From  (97),  we  outline  the  equations  for  Coio  and  C0n 


dC'oio 

dt 


+  hCon  —  2eCooo  —  A^Guo  +  E  A(E,iG-iio? 

rl^O 


dCpu 

dt 


-  hCoio  =  A ooiCon  +  E  Ajri,iCrin. 

r  1^0 


(101) 


Here,  Cooo  does  not  depend  on  time  and  set  to  be  equal  unity.  In  a  system  (101) 
under  sign  of  summation,  the  moments  Cr \m  enter  at  m  —  0, 1  and  r  /  0,  and,  in 
general  case,  to  calculate  the  temporal  dependence  of  a  drift  velocity,  a  full  system 
of  the  moment  equations  is  needed  to  be  solved. 

In  the  case  of  the  diagonal  interaction  matrix  ,  the  equations  for  Coio  and  C011 
are  split  out  and,  taking  in  mind  (98)  from  (101),  we  have 


+  kCQ 


2e  —  —On 


—  hCn  in  —  —  On 


(102) 


Let  till  t  —  0  there  was  a  magnetic  field  h  and  no  electric  field,  and  at  t  —  0  the 
electric  field  e  experiences  a  sharp  switch-on. 

As  the  magnetic  field  itself  does  not  affect  the  equilibrium  DF,  the  DF  continues 
to  be  Maxwellian  till  a  switch-on  of  the  electric  field/,  and  C'0;m  =  0  at  l  ^  0.  Hence,  a 
system  (102)  has  zero  initial  conditions.  Apply  the  Laplace  perform  gp  =  L(g(t))  = 

OO 

/  e~ptg{t)dt  to  (102) 
o 

ip  +  1)^010  +  hCp0U  =  — 

P 

— hCoio  +  ip  +  l)Coii  =  0.  (103) 

This  system  has  a  determinant  D  =  (p  +  l)2  +  h2  and  its  solution  takes  a  form 


np  2(p+l)e  ^  2  he 

G°10~  pD  ’  C°"-pD- 

Carrying  out  the  inverse  Laplace  perform,  we  obtain 


(104) 


C'oio(i)  =  ~  +^9  (l  +  e  *  ih  sin  ht  -  cos  ht)j , 


Con(^)  — 


1  +  h2 


sin  ht 


+  cos  ht 


(105) 
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Figure  7:  The  temporal  dependences  of  the  moments  Coio  and  Con  for  arbitrary 
Maxwellian  molecules,  e  =  1,  h  —  1  in  the  non-standard  normalization.  Solid  curves 
present  the  analytical  results,  points  are  the  numerical  solutions  of  the  system  of 
the  moment  equation.  Figure  a  —  Cn0,  b  —  Cm 


The  comparison  of  these  analytical  results  with  those  obtained  via  the  numerical 
solution  of  the  system  of  the  moment  equations  at  e  =  1  and  h  =  1  is  presented 
in  Fig.  7.  One  can  see  that  the  temporal  dependences  have  an  oscillating  compo¬ 
nent.  When  h  increases,  the  oscillation  feature  of  the  transient  process  evolves  more 
evidently. 

At  the  stationary  state,  from  (105)  at  t  — >  oo  we  obtain 


2e 

“  1  +  K2' 


Con  — 


2  he 

1  +  h2' 


(106) 


In  papers  [33],  [8],  an  expression  for  a  drift  velocity  at  the  stationary  state  is  con¬ 
sidered  for  the  constant  collision  frequency,  i.e.,  at  the  diagonal  interaction  matrix. 

If  we  transfer  to  the  our  notations,  the  results  [33]  are  totally  coincide  with 
(106).  At  the  stationary  state,  the  Lorentz  angle  (99)  and  the  drift  velocity  (100)are 
determined  as  follows 

tan  ip  =  h, 
e  e 

ud  = - — r  =  ecos^  =  -sinip.  (107) 

(1  +  h2)  2  ii 

If  we  go  to  the  dimensional  variables,  it  is  seen  that  the  oscillation  frequency  of 
the  drift  velocity  coincides  with  the  Larmor  one  in  the  non-stationary  process. 


In  the  7th  quarterly  report,  the  expression  of  the  principal  physical  moments  is 
obtained  via  the  moments  Crim  and  it  is  shown  that  for  their  calculation  Cr/m  with 
l  <  2  and  r  <  1  with  are  needed.  These  formulas  were  obtained  in  the  quarterly 
report  and  do  not  presented  here.  Besides  a  drift  velocity,  there  are  considered  a 
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pressure  tensor  and  a  vector  of  a  heat  flow  as  the  physical  moments;  a  pressure 
tensor  has  a  component  pxz  besides  of  the  diagonal  components  pxx ,  pyy  and  pzz. 
In  the  case  of  the  diagonal  interaction  matrix,  further  splitting  out  occurs  of  the 
system  of  the  moment  equations  for  Crim.  At  fixed  r  and  l,  there  are  l  + 1  equations 
for  0  <  m  <  l.  These  equations  are  solved,  also,  using  the  Laplace  transform  and 
convolution  theorem. 


Figure  8:  The  trajectories  of  a  drift  velocity  for  several  values  of  h.  Curve  1  refers 
to  h=. 5,  2  -  h  =  1,  3  -  h  =  2,  4  -  h  =  5,  5  -  h  =  10. 

As  a  result,  the  temporal  dependences  of  all  physical  moments  are  built  up.  The 
mobility  vector  behaviour  is  of  interest.  As  one  can  see  in  Fig.  8,  it  oscillates  and, 
on  a  (K,  Khai)- plane,  it  represents  a  twisting  spiral,  the  oscillation  amplitude  and 
the  number  of  the  spiral  turns  increase. 

Thus,  for  the  first  time,  Report  [32]  presents  the  investigations  of  the  non-  sta¬ 
tionary  transient  processes  evolving  into  an  ion  admixture  on  the  background  of  the 
atoms  at  a  sharp  switch-on  of  the  electric  held  under  the  magnetic  held  (generally, 
for  the  pseudo-Maxwellian  molecules). 

A  code  was  created  to  solve  a  system  of  the  integro-differential  equations  involv¬ 
ing  the  ion  interaction  with  the  electric  held  as  well  as  the  magnetic  one.  It  was 
shown  that  the  system  of  the  moment  equations  is  divided  on  the  subsystems  for 
Crim  with  hxed  r  and  l  and  0  <  m  <  l. 

Involving  the  Laplace  transform,  this  subsystem  of  the  differential  equations  is 
converted  into  a  subsystem  of  algebraic  equations.  It  was  proved  that  the  determi¬ 
nant  of  this  subsystem  has  the  roots  A ri  ±  ikh ,  where  k  are  integers  of  which  parity 
coincides  with  that  of  l. 

It  was  maintained  that,  in  the  case  of  the  diagonal  interaction  matrix  (Maxwellian 
molecules  with  arbitrary  angular  part  of  the  scattering  cross  section),  the  subsys¬ 
tems  of  the  equation  aligned  in  a  recurrently  solvable  sequence.  It  signifies  that,  if 
one  moves  along  this  sequence,  solving  the  subsystem  for  Cxlm,  all  Cpr,Urn,  with  r'V 
combinations  distinct  from  r’V  were  already  known. 
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The  analytical  solutions  for  the  moments  Cr[m  with  the  small  values  of  r,  /.  It 
was  shown  in  what  manner  were  expressed  the  principal  physical  moments  via  these 
moments.  It  was  found  that,  the  oscillations  arose  during  the  transient  stage  of  the 
process  evolution.  The  principal  harmonic  at  a  Larmor  frequency  as  well  as  higher 
harmonics  at  the  frequencies,  which  are  multiples  of  the  principal  one,  were  excited. 

Of  exceptional  interest  is  excitation  of  oscillations  of  the  Hall  current  and  a 
current  along  the  electric  field  which  can  be  observed  in  a  particular  experiment 
and  use  for  diagnostics. 

Via  the  analytical  formulas  as  well  as  the  numerical  solution  of  the  system  of 
the  moment  equations  with  aid  of  our  code,  the  large  amount  of  the  results  were 
obtained  to  describe  evolution  of  both  the  DF  moments  evolution  and  physical  ones. 
The  fine  coincidence  of  the  analytical  and  numerical  results  was  obtained. 

A  detailed  analysis  of  a  behaviour  of  physical  moments  is  presented  at  the  tran¬ 
sient  stage  under  different  values  of  the  magnetic  field  h  and  dependences  of  the 
stationary  values  of  these  moments  on  h. 


8  Eighth  quarter 


The  8th  quarterly  report  [34]  describes  the  calculations  of  the  ion  physical  moments 
at  the  mutually  perpendicular  electric  and  magnetic  fields  for  various  interaction 
models.  As  before,  under  investigations  is  a  transient  process  after  the  instant 
switch-on  of  the  electric  field.  The  7th  quarterly  report  [32]  develops  the  code  for 
solving  the  system  of  the  moment  equations  in  the  case  arbitrary  interaction  model. 
It  is  shown  that  the  numerical  and  analytical  results  coincide  in  practice,  for  the 
Maxwellian  molecules,  this  proving  the  involved  code. 

The  8th  quarterly  report  uses  widely  the  developed  code  to  study  the  transient 
processes  at  the  various  ion-atom  interaction  models.  When  calculating,  the  non¬ 
standard  normalization  was  used. 

First  of  all,  the  calculations  for  the  interaction  models  with  the  constant  collision 
frequency  (in  this  case,  the  interaction  matrix  is  diagonal)  and  a  comparison  of  the 
processes  for  the  pseudo-Maxwellian  molecules  and  CEM-model  is  undertaken.  It 
is  shown  that,  in  the  case  of  the  diagonal  interaction  matrix,  the  physical  moments 
are  proportional  to  the  different  degrees  of  the  dimensionless  electric  field  e. 


K=^L, 


A  hril 


^ dx  —  Pxx  A  _  Pyy  A  _  Pzz  A 
)  Pxx  =  2  ’  Pyy  =  2  ’  Pzz  =  2  ’ 


Pxz  =  *f,  q  =  i,  (108) 

ez  e 6 

which  do  not  depend,  already,  on  e. 

For  all  models  with  the  non-diagonal  interaction  model  in  a  domain  of  the  weak 
electric  field,  the  physical  moments  turn  out  to  be  proportional  the  different  degrees 
of  e.  However,  if  a  drift  velocity  (mobility)  and  pressure  tensor  depend  on  e  in  the 
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manner  as  in  the  case  of  the  diagonal  matrix,  a  heat  flow  turns  out  to  be  proportional 
to  e  but  not  e3. 

Besides  the  models  with  r  =  const,  the  systematic  calculations  was  carried 
out  for  the  following  interaction  models:  1.  HS-model,  the  hard-sphere  model.  2. 
CEHS-model,  the  charge  exchange  with  the  scattering  cross  section  independent  on 
the  velocity.  3.  Coulomb  interaction.  For  all  models,  the  reduced  physical  moments 
are  determined  as  follows 


K  = 


Udz 


h  hal 


H dx  —  Pxx  -5  _  Pyy  —  Pzz  -h 

— — ;  Pxx  =  “5  >  Pyy  =  ~2  ’  Pzz  =  72  ; 


Pxz 


Px 


-  q 

q=  - 

e 


(109) 


We  call  the  weak  electric  held  that  of  e  <  0.1.  In  a  subregion  of  the  weak 
held  e  <  10~2,  a  change  in  e,  for  all  models  under  consideration  does  not  results  in 
any  changes  in  the  temporal  dependences  of  the  reduced  physical  moments.  The 
reduced  physical  moments  at  the  very  weak  electric  hels  (e  <  10~2)  are  called  by  us 
the  standard  reduced  physical  moments. 

For  all  model  except  r  =  const,  with  an  increase  in  the  electric  held,  the  reduced 
physical  moments  experience  the  changes.  Nevertheless,  the  reduced  physical  mo¬ 
ments  turn  out  to  be  advisable  to  study  when  being  built  up  not  only  at  the  small 
e.  The  8th  quarterly  report  gives  the  results  for  three  interaction  models,  5  values 
of  e  variants,  15  overall,  being  considered  in  the  domains  of  the  weak  and  moderate 
electric  held.  In  every  variant,  the  temporal  dependences  of  9  physical  moments  are 
presented,  being  built  up  at  6  values  of  h  (h  =  0.1,  h  =  0.5,  h  =  1,  h  =  2,  h  =  5, 
h  =  10).  The  total  amount  of  dependences  is  above  800.  All  this  set  of  the  curves 
can  be  called  ATLAS  of  the  temporal  dependences  of  the  physical  moments  for 
the  ions  at  the  crossed  electric  and  magnetic  hclds  under  a  sharp  switch-on  of  the 
electric  held. 

This  atlas  is  the  principal  part  of  the  full  report  for  the  Quarter  8. 

It  was  found  that,  at  the  weak  magnetic  held  in  the  domains  of  the  moderate  and 
strong  electric  held,  a  distinction  from  the  reduced  physical  moments  is  pronounced, 
being  especially  large  for  the  heat  how.  However,  with  an  increase  in  the  magnetic 
held,  this  distinction  diminishes  and,  at  h  —  5,  10,  the  reduced  physical  moments 
coincide  with  the  standard  ones.  At  the  weak  magnetic  held  (h  =  0.1,  h  =  0.5),  for 
the  models  A  =  const  both  the  mobility  and  the  tensor  components  turn  out  to  be 
smaller  that  the  standard  ones.  A  completely  different  situation  occurs  in  the  case 
of  the  reduced  physical  moments  in  the  case  of  the  Coulomb  interaction.  In  Fig.  9, 
the  mobility  behaviour  is  shown. 

For  this  model,  the  mobilities  turn  out  to  be  higher  than  the  standard  ones. 
Moreover,  at  a  small  h  and  e  =  0.5,  there  is  a  tendency  to  an  unlimited  increase;  at 
the  smaller  values  of  e,  this  effect  does  not  almost  occur,  but,  at  e  >  0.5,  it  is  even 
more  pronounced.  Along  an  analogy  with  the  electron  running-off,  we  call  it  as  ion 
running-off. 
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Figure  9:  Dependences  K(t)  (Curve  1)  and  Khai(t )  (Curve  2)  for  the  Coulomb 
interaction  at  e  =  0.5  and  different  values  of  h.  The  dashed  curves  refer  to  the 
standard  mobility.  Figure  1  a  -  h  =  .1,  b  -h  =  .5,  c  -h  —  1,  d  -  h  =  2,  e  -h  —  5,  / 
-h  =  10. 

This  effect  is  a  consequence  of  a  very  strong  decrease  in  the  scattering  cross 
section  with  an  increase  in  a  relative  velocity(cr  ~  g~A)- 

From  Fig.  9,  one  can  see,  also,  that  the  running-off  effect  is  suppressed  by  the 
magnetic  field.  At  the  strong  magnetic  held  (h  =  5  or  h  —  10)  the  reduced  mobility 
coincides  with  the  standard  one,  the  moment  method  convergence  being  restored.  A 
good  convergence  at  the  strong  magnetic  held  is  observed  for  all  interaction  models, 
it  makes  possible,  at  large  h,  an  advance  into  the  strong  electric  held  domain,  i.e.,  in 
the  very  region  being  non-calculable  with  standard  moment  method  at  no  magnetic 
held. 

Thus,  in  Quarter  8,  the  calculation  of  the  evolution  process  concerning  the  phys¬ 
ical  moments  for  various  laws  of  interaction  from  the  most  hard  interaction  potential 
for  the  hard-sphere  model  to  the  softest,  i.e.,  Coulomb  interaction. 

Below,  the  conclusions  made  on  a  base  of  this  study: 

1.  The  notion  of  the  reduced  physical  moments  (see,  formulas  (108)  and  (109)) 
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is  introduced.  The  reduced  moments  in  a  domain  of  the  weak  field  do  not  depend 
on  e  and  are  used  as  the  standard  reduced  moments. 

2.  It  is  found  that  the  behaviour  of  the  components  of  a  pressure  tensor  is 
determined,  generally,  by  the  angular  part  of  the  scattering  cross  section  and  that 
of  the  components  of  a  heat  flow  by  the  velocity  dependence  of  the  scattering  cross 
section. 

3.  It  is  shown  that,  in  a  domain  of  the  weak  electric  field,  the  mobilities  in 
a  direction  of  the  electric  field  K  as  well  as  that  perpendicular  to  the  electric  and 
magnetic  field,  Khai,  are  similar  in  their  behaviour  at  all  the  scattering  cross  sections. 

4.  In  a  domain  of  the  moderate  and  strong  electric  fields,  for  a  series  of  the 
interaction  cross  section,  the  ATLAS  of  the  non-stationary  characteristics  of  the 
physical  moments  is  compiled,  these  characteristics  are  presented  in  a  wide  range  of 
e  and  h  changes.  5.  It  is  shown  that,  in  the  case  of  the  Coulomb  interaction  at  a 
small  h,  the  ion  running-off  takes  place  and  this  process  becomes  more  pronounced 
with  an  increase  in  e. 

The  running-off  process  at  the  Coulomb  interaction  is  suppressed  by  the  strong 
magnetic  field. 

6.  It  is  shown  that,  for  any  interaction  models  at  the  strong  magnetic  field,  all 
the  physical  moments  coincide  with  the  standard  ones. 

7.  The  moment  method  can  be  used  successfully  for  calculation  of  the  non- 
stationary  ion  process  at  the  crossed  electric  and  magnetic  fields  in  a  domain  of 
the  very  strong  electric  fields  if  the  magnetic  field  is  sufficiently  strong(e  and  h  are 
commensurable) . 


9  Ninth  quarter 

In  the  9th  quarterly  report  [35],  the  calculation  of  the  ion  DF  evolution  at  the 
mutually  perpendicular  electric  and  magnetic  fields  were  performed  for  different 
interaction  models.  The  calculations  were  based  on  a  code  concerned  with  evolution 
of  the  moments  Crim  which  was  developed  in  Quarter  7  [32],  The  problems  solved 
in  Quarters  9  and  10  involve  the  DF  restoring  via  its  moments  as  well  as  descriptive 
DF  representation. 

With  no  magnetic  field,  the  electric  field  switch-on  directed  along  z  results  in 
stretching  the  DF  along  cz.  The  DF,  at  any  time,  has  the  axial  sysmmetry  relative 
to  c^-axis,  i.e.,  depends  on  two  variables  cz-  and  cp  =  ^/c|  +  cy. 

There  are  more  complication  when  dealing  with  the  DF  shape  as  well  as  the 
descriptive  representation  of  its  evolution  at  the  switch-on  of  the  mutually  perpen¬ 
dicular  electric  and  magnetic  field.  A  direction  of  the  magnetic  field  is  matched  with 
y-axis,  Due  to  ion  rotation  about  y  axis,  there  is  a  strong  rebuilding  of  the  DF  in 
(cXjCz)- plane,  the  symmetry  relative  cy  substitution  with  —  cy  being  conserved.  A 
departure  from  the  maxwcllian  dependence  relative  to  \cy\  turns  out  to  be  neglecting. 
So,  we  follow  the  DF  evolution  when  building  up  f(cz,  cx)  at  cy  =  0. 
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At  such  image,  in  a  fixed  time,  the  DF  is  represented  by  certain  surface  in  3D 
space.  The  DF  expansion  can  be  performed  in  two  steps:  first,  an  expansion  in 
terms  of  the  spherical  harmonics: 

/(c)  =  ^FZm(0,^)/Zm(c),  (110) 

l,m 

second,  in  terms  of  the  Sonine  polynomials 

OO 

fim(c )  =  M(T0,  c2)  £  CV,*>mc^S/2(c2)  (m) 

r=0 

In  the  problems  under  consideration,  the  expansion  is  performed  in  terms  of  the 
cosine  spherical  harmonics,  so,  in  (110),  an  index  i  is  omitted.  Knowing  Crim,  using 
(111),  one  can  build  up  all  fim(c),  then,  using  (110),  restore  the  DF  /(c). 

If  e  <  1,  the  DF  can  be  expanded  in  a  series  in  terms  of  the  powers  of  e.  As  a 
preliminary,  we  extract  the  Maxwellian  distribution  from  the  DF,  passing  to  /(c) 
via  a  formula 

/ (c)  =  M (T0,  c2)  +  / (c)  (112) 

Clearly,  for  the  expansion  coefficients  of  the  spherical  harmonics  we  have 

Zoo(c)  =  /oo(c)  -  M(T0,  c2),  if  l  =  0,  /{m(c)  =  /{m(c)  if  l  ^  0  (113) 

In  the  9th  quarterly  report,  the  dependences  fim  on  e  were  studied  numerically  in 
a  region  e  <  1  at  the  mutually  perpendicular  electric  and  magnetic  fields,  It  was 
found  that,  independently  on  the  interaction  model,  fim  ~  el  at  l  0  and  f00  ~  e2. 
The  functions  <f/m  were  determined  as  follows 

^  /  >  0, 

e 

$oo  =  %  (114) 

In  a  domain  of  the  very  weak  electric  field,  these  functions  do  not  depend  on  e. 
The  calculations  of  <3>;m  functions  were  carried  out  for  different  interaction  mod¬ 
els.  Graphically,  they  do  not  differ  at  e  <  10~2.  Under  an  increase  in  e,  these 
function  change.  The  report  presents  the  evolution  of  several  such  functions  at 
different  values  of  the  electric  and  magnetic  field. 

In  the  case  of  the  weak  electric  field,  when  expanding  the  DF  in  terms  of  e 
powers,  there  can  be  a  restriction  by  the  first  term 

/(c)=M(T0,c2)(l  +  e$(c))  (115) 

A  function  <3>(c)  independent  on  e  in  a  domain  of  the  weak  electric  field  is  called  the 
universal  DF,  being  under  study  in  the  report.  Fot  it,  we  have 

i 

$(c)  =  £  $1m(c)Fim(0,  (p),  (116) 

m= 0 
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A  polar  angle  p  in  (cx,  cy)-plane  is  counted  from  the  direction  of  cx-axis.  Hence,  in 
cy  =  0-plane,  an  angle  p  =  0  and 

$(c*,  Cz)  =  $io(c)  cos(0)  +  <hn (c)  sin(0).  (117) 

Denote 

m  =  p<S.l0(c)?  +  (tn(c)y.  (118) 

and  determine  an  angle  a(c)  as  follows 

(c) 

a(c )  =  arctan  — — —  (119) 

^io(c) 

Then,  $n(c)  =  R(c )  sina(c)  and 

®(cx,  cz)  =  R{c)  ■  cos(0  -  a(c)).  (120) 

From  (120)  one  xan  see  that  the  DF  $(cx,c~)  is  uniquely  determined  by  two 
functions  R(c)  and  a(c)  which  are  called  its  characteristics.  At  hxed  c,  the  maximum 
of  $(0)  is  observed  at  0  =  a,  being  equal  R(c),  and  the  minimum  being  at  0  = 
a  +  7T.  Additionally,  we  have 


®(-cx,  -cz)  =  -$(cx,  cz) 


(121) 


Figure  10:  The  DF  $(cx,c2)  at  e  <  10-2  and  its  characteristics,  CEM-model,  h  =  0. 
Figure  a  is  3D-  representation  of  the  stationary  DF  (t  — >  oo).  Figure  b  represents 
the  characteristics  R(c)  at  different  times,  Curve  1  refers  to  t  —  0.1,  2  -  t  —  0.5,  3  - 
t  —  1,  4  -  t  —  3,  5--—  t  =  5. 
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a) 


Figure  11:  The  DF  $(cx,cz)  at  e  <  10~2  and  it’s  characteristics,  CEM-  model, 
h  =  1.  Figure  a  is  3D-  representation  of  the  stationary  DF  (t  — >  oo).  Figure  b  and 
Figure  c  are  projection  of  a-linc  at  the  cx,cz- plane.  Curve  1  refers  to  t  =  0.1,  2  - 
t  =  0.5,  3  -  t  =  1,  4  -  t  =  3,  5  -  t  =  5. 


Figure  12:  The  DF  <&(cx,  cz)  at  e  <  10-2  and  it’s  characteristics,  HS-model,  h  =  0. 
Figure  a  refers  to  3D-representation  of  the  stationary  DF  (t  — >■  oo).  Figure  b  presents 
the  characteristics  R(c)  at  different  times.  Curve  1  refers  to  t  —  0.1,  2  -  t  —  0.5,  3 
-t  =  l,  4-t  =  3,  5-t  =  5. 

This  property  gives  a  possibility  to  be  restricted  by  the  DF  building-up  only  at 
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a  half-plane,  e.g.,  cz  >  0. 

A  dependence  a(c)  is  a  projection  of  a  line  of  maxima  of  the  DF  in  the  polar 
coordinates  on  (cx,  c2)-plane. 

It  should  be  noted  that,  with  no  magnetic  field,  a(c)  =0. 

When  graphically  representing  the  DF,  3D  picture  is  built  up  in  which  there  is 
&(cx,Cz)  and  a  thick  line  represents  the  line  of  maxima.  3D  figure  is  built  up  for 
the  stationary  state  corresponding  to  t  — >  oo.  For  obtaining  a  representation  of 
evolution  process,  additionally,  for  several  times,  the  characteristics  of  a  function 
d>(cx,  cz)  are  shown,  i.e.,  the  curves  R(c )  and  projections  a- line  at  a  plane  (cx,cz); 
these  projections  are  not  built  up  at  no  magnetic  field  (h  —  0). 

Several  examples  of  the  calculations  are  shown  in  Figs.  10  -  13 


Figure  13:  The  DF  fI>(cx,  cz)  at  e  <  10~2  and  its  characteristics,  HS-model,  h  =  1. 
Figure  a  refers  to  3D-representation  of  the  stationary  DF  (t  — >  oo).  Figures  b  and 
c  present  the  DF  characteristics  at  different  times.  Figure  b  refers  to  R(c),  Figure 
c  presents  the  projections  of  a-line  at  a  palne  cx,cz.  Curve  1  refers  to  t  —  0.1,  2  - 
t  =  0.5,  3  -  t  =  1,  4  -  t  =  3,  5  -  t  =  5. 


43 


10  Tenth  quarter 

In  the  10th  quarterly  report  [36],  the  calculations  of  DF  evolution  were  continued 
for  different  interaction  models  at  the  instant  switch-on  of  the  electric  field  and  the 
mutually  perpendicular  electric  and  magnetic  fields,  Being  distinct  from  Quarter  9 
[35],  when  the  electric  field  was  set  to  be  weak  (e  <  10”2),  here,  it  was  considered 
the  significantly  stronger  electric  field  up  to  e  =  1  and  even  more.  As  in  the  previous 
report,  the  expansion  in  terms  of  the  real  spherical  harmonics  was  of  wide  use. 

If,  under  the  weak  electric  field  at  fixed  velocity  c,  the  DF  dependence  on  an 
angle  of  rotation  in  (cx,  cQ-plane  is  of  cos-character,  now,  this  dependence  is  more 
complicated  and  changes  when  changing  both  c  and  magnetic  field  h.  It  changes 
strongly  when  changing  the  interaction  model. 

An  increase  in  the  electric  field  with  no  magnetic  field  results  in  stretching  the 
DF  along  c2-axis,  and  a  sharp  maximum  arises  on  the  angular  dependence.  At  the 
no  so  strong  magnetic  field,  the  sharp  maximum  is  conserved,  a  change  is  only  a 
direction  of  the  stretching.  The  maximum  height  can  experience  a  change  of  several 
orders  of  value  when  changing  c.  In  this  case,  all  that  results  in  a  loss  of  clearness  of 
3D-representation  of  the  DF,  which  was  used  in  the  9th  quarterly  report  [35]  when 
studying  the  DF  at  small  e. 

So,  to  represent  the  DF,  the  angular  dependences  in  (cx,cz)~ plane  were  built  up 
at  several  values  of  a  velocity  module  c  =  +  c2.  An  angle  i/j  was  counted  from 

c2-axis  and  changed  from  — pi  to  pi. 

In  the  10th  quarterly  report,  the  detailed  calculation  were  carried  out,  being 
illustrated  by  a  number  of  the  figures  for  three  interaction  models:  HS-model,  CEM- 
model,  and  the  Coulomb  interaction.  Here,  as  an  example,  we  present  only  four 
figures  (  14-  17),  corresponding  to  two  extreme  cases:  the  hardest  model  (HS- 
model)  and  the  softest  (Coulomb  interaction).  Only  data  are  presented  for  one 
value  of  the  electric  field,  e  =  1,  and  two  values  of  the  magnetic  field,  h  —  0  and 
h  =  1. 

In  these  figures,  one  can  see  that  the  magnetic  held  results  in  a  rotation  of  the 
location  of  the  maximum,  being  more  pronounced  in  a  domain  of  the  high  velocities. 
Thus,  e.g.,  a  transfer  from  HS-model  to  the  Coulomb  interaction  results  in  about 
two  order  higher  maximum  at  c  =  5. 

In  the  case  of  the  Coulomb  interaction,  a  point  under  consideration  was  in  what 
manner  the  running-off  influences  the  DF  evolution.  This  effect  is  observed  at  not  so 
strong  transverse  magnetic  held  and  appears  in  a  constant  and  rather  rapid  increase 
of  the  maximum  of  a  dependence  ffy)  in  a  domain  of  high  velocities.  This  problem 
is  rather  complicated  and,  when  using  the  standard  moment  method,  we  come  up 
quite  quickly  (at  not  so  long  time)  against  a  situation  of  no  convergence  of  the  DF. 

To  ameliorate  the  convergence,  a  transfer  to  other  basis  was  undertaken,  i.e., 
instead  of  an  expansion  about  the  weight  Maxwellian  with  a  temperature  of  atoms 
it  was  used  the  Maxwellian  of  higher  temperature. 
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Figure  14:  The  dependences  at  different  times  at  different  values  of  the  velocity 
module  c  for  HS-model;  e  =  1,  h  =  0.  Figure  a)  refers  to  c  =  0.5;  b)  -  c  =  1;  c)  - 
c  =  3;  d)  -  c  =  5. 

Represent  the  DF  in  two  bases  of  the  temperatures  T0  and  Tf: 

OO  OO 

/(»,  t)  =  EE  C°j(t)M(v,  T0)S^(v2/T0)c‘P,  = 

1=0  k=0 

oo  oo 

=  E  £  0)V(”,  Tl)S[n{v2IT1)<iP,.  (122) 

1=0  r= 0 

The  vectors  C°  and  C1,  when  transferring  from  a  basis  T0  to  a  basis  Tj  ,  change  in 
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a)  f  b) 


Figure  15:  The  dependences  f{ip)  for  different  times  at  different  values  of  the  velocity 
module  c  for  HS-model;  e  =  1,  h  =  1.  Figure  a)  refers  to  c  =  0.5;  b)  -  c  =  1;  c)  - 
c  =  3;  d)  -  c  =  5. 


a  manner  as  follows 

OO 

C\i  =  Y,Dr,k,i{TuTo)C\ 

k= o 

A  transfer  matrix  Dr^j  is  as  follows 


o 


Dr,k,i  — 


(-1  )r+kk\  To+l/2  (Ti  -  T0)k~ 

y\  rjik+l/2  ^ 


(123) 


(124) 


Denote  a  ratio  of  the  temperature  of  these  bases  as  RT  =  T\  /T0.  We  call  the 
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a)  f  b) 


Figure  16:  The  dependences  f{ip)  for  different  times  at  different  values  of  the  velocity 
module  c  for  the  Coulomb  interaction;  e  =  1,  h  =  0.  Figure  a)  refers  to  c  =  0.5;  b) 
-  c  =  1;  c)  -  c  =  3;  d)  -  c  =  5. 

modified  moment  method  that  in  which  a  parameter  RT  changes  in  time.  A  point 
concerning  the  optimum  choice  of  a  temporal  dependence  of  a  parameter  RT  at 
different  e  requires  an  additional  study. 

A  simpler  method  was  used  called  ” semi- modified  moment  method”.  In  this 
method,  the  calculations  of  the  moment  are  carried  out  in  an  initial  basis  (temper¬ 
ature  of  atoms)  and,  before  the  DF  building-up,  a  transfer  to  a  new  basis  is  caaried 
out. 

It  occurs  that,  using  this  method  we  advance  further  in  studying  the  running-off 
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a)  f  b) 


c)  f  d) 


Figure  17:  The  dependences  f{ip)  at  different  times  at  different  values  of  the  velocity 
module  c  for  Coulomb  interaction;  e  —  1,  h  —  1.  Figure  a)  refers  to  c  =  0.5;  b)  - 
c  =  1;  c)  -  c  =  3;  d)  -  c  =  5. 

effect  in  the  case  of  the  Coulomb  interaction.  At  c  =  5,  the  transfer  to  a  new  basis 
being  over,  we  succeed  in  description  of  the  running-off  process  with  about  one  order 
increase  of  the  maximum  on  the  angular  dependence  of  the  DF. 

Thus,  in  Quarter  10,  the  DF  calculations  for  three  interaction  models  were  per¬ 
formed  at  the  different  magnetic  fields  in  the  domains  of  the  moderate  and  strong 
electric  field.  Onr  investigations  being  performed,  the  results  were  obtained  as  fol¬ 
lows 

1.  The  DF  evolution  is  strongly  dependent  on  the  features  of  particle  interaction. 
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Especially  strongest  is  an  influence  of  interaction  potential  on  the  DF  in  a  domain 
of  high  velocities. 

2.  At  zero  magnetic  field,  in  a  domain  of  high  velocities,  in  the  vicinity  -0  =  0 
there  is  maximum  of  which  height,  at  the  same  electric  field,  becomes  enormously 
higher  when  there  is  a  transfer  of  the  hard  potential  (HS-model)  to  the  softer  one 
(CEM-model). 

3.  With  an  increase  in  the  magnetic  field,  there  is  a  rotation  of  location  of  the 
maximum  toward  positive  0,  and  the  maximum  height  decreases;  in  a  domain  of 
the  strong  magnetic  fields  h  >  1,  relaxation  process  of  the  maximum  becomes  non- 
monotonous. 

4.  At  the  very  strong  magnetic  field  (h  =  5,  h  =  10),  the  DF  evolution  does  not 
depend,  in  practice,  on  a  mode  of  the  particle  interaction  but  it  is  determined  by 
the  magnetic  field. 

5.  At  the  DF  evolution  in  the  case  of  the  Coulomb  interaction  and  the  small 
values  of  h,  a  general  process  is  the  ion  running-off.  In  a  domain  of  large  c,  the  DF 
increases  monotonously  with  no  saturation. 

6.  To  study  the  ion  running-off,  the  semi-modified  moment  method  was  de¬ 
veloped  with  the  calculation  of  the  moments  Cri  in  an  initial  basis,  and,  before 
calculation  of  the  DF,  there  is  a  transfer  to  a  new  basis  with  higher  temperature. 

7.  In  the  case  of  the  Coulomb  interaction  and  the  very  strong  magnetic  field 
in  a  domain  of  high  velocities,  unusual  phenomena  were  revealed  related  with  the 
intense  oscillations  of  a  cloud  of  the  fast  particles. 


11  Eleventh  quarter 

The  11th  quarterly  report  [37]  deals  with  ion  diffusion  at  the  electric  field. 

A  small  ion  admixture  being  presented  at  the  background  of  a  number  of  atoms, 
only  transport  process  exists,  i.e.,  diffusion,  besides  the  transfer  of  charges  under  the 
electric  field.  To  study  this  process,  we  use,  also,  the  non-  stationary  approach.  With 
this  aim,  at  t  —  0,  a  concentration  gradient  Vn  is  artificially  introduced  into  initially 
uniform  ion  ensemble.  This  gradient  is  supposed  to  be  constant  in  a  neighborhood 
of  a  point  r0  under  consideration,  the  size  of  this  neighbourhood  should  be  many 
times  the  free  length  path.  Solving  the  Boltzmann  equation,  the  DF  evolution  is 
found  as  well  as  diffusion  coefficient.  At  sufficiently  long  time,  the  process  saturation 
occurs  and  a  stationary  value  of  the  diffusion  coefficient  is  found  and  corresponding 
DF.  The  solution  of  the  Boltzmann  equation  is  built  up  via  the  moment  method 
with  a  number  of  the  moments  as  well  as  use  of  all  store  of  knowledge  on  the  ion 
mobility  calculation.  At  the  weak  electric  field,  mobility  and  diffusion  effects  do 
not  depend  of  each  other.  At  the  more  strong  electric  field,  the  diffusion  coefficient 
depends  on  the  electric  field  value.  In  this  case,  two  diffusion  coefficients  are  built  in: 
transverse  and  longitudinal.  Independently  on  the  electric  field  value,  as  it  usually 
proceeds  [4],  the  concentration  gradient  value  is  set  as  the  small  one.  In  this  case, 


49 


only  the  first  order  perturbation  theory  can  be  used.  Under  the  sufficiently  high 
electric  field,  the  problem  is  solved  as  follows.  First,  using  the  developed  codes,  the 
problem  concerning  the  ion  evolution  after  a  sharp  switch-on  of  the  electric  field  of 
given  value  is  solved.  The  stationary  DF  obtained  in  this  problem  is  chosen  as  an 
initial  condition  to  solve  the  next  problem  concerning  the  DF  evolution  after  the 
sharp  switch-on  of  Vn. 

In  the  11th  quarterly  report,  the  known  results  concerning  the  ion  diffusion  and 
mobility  are  described  in  details  as  well  as  the  various  variants  of  the  transfer  to  the 
dimensionless  variables.  It  is  shown  that  a  dimensionless  electric  held  €m  from  [4] 
is  four  times  e  determined  by  us.  Hence,  6m  =  4  corresponds  to  (e  =  1)  the  latter 
being  the  low  boundary  of  the  strong  held  determined  by  us. 

In  the  dimensionless  formulation,  the  Boltzmann  equation  for  ion  diffusion  at  no 
electric  held  takes  a  form 

^  +  6cJ  =  j(f,M ),  (125) 

and,  for  ion  motion  at  the  electric  held  (with  no  diffusion  due  th  concentration 
gradient)  it  is  as  follows 

f  +  ef  =  U.M).  (126) 

where  J(f,M )  is  the  dimensionless  collision  integral, 


dn/dq  dn/dz 

o  = - =  A - , 

n  n 


(127) 


Here,  a  dimensional  coordinate  q  =  z/X  is  used. 

Let,  in  the  problem  concerning  the  mobility,  we  have  the  very  weak  electric  held 
(e  -C  1),  and  the  very  small  density  gradient  (5  <C  1)  is  involved  in  the  problem 
concerning  diffusion.  In  Report  [28],  the  case  of  the  weak  electric  held  is  treated 
in  details.  If  (5  -C  1)  and  (e  -C  1),  the  DF  is  sought  in  a  form  fs(c )  =  M(c)(l  + 
5(fs( c)),  referring  to  the  Chapman- Enskog  method,  and  fe( c)  =  M(c)(l  +  etpe(c)), 
correspondingly. 

If  the  terms  of  the  second  and  higher  orders  are  neglected  and,  in  the  problem 
concerning  mobility,  we  write  dM/dcz  =  —2 czM  for  the  derivative  with  respect  to 
the  Maxwellian,  from  (125)  and  (126),  it  follows 


M~^t  +  °zM  = 


(128) 


M -  2 czM  =  J(Mipe ,  M). 
Comparing  (129)  and  (128),  we  conclude  that 


(129) 


<Ps(c,t)  =  -~<pe(c,t) 


(130) 
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It  follows  that,  for  the  mean  dimensionless  velocities  in  the  cases  under  consideration 
we  have 

1 

US  =  --Ue.  (131) 

The  dimensionless  mobility  K  and  diffusion  coefficient  D  are  determined  as  follows 

ue  =  Ke,  us  =  — D5 .  (132) 

From  (131)  and  (132),  it  follows 

D  =  i  K.  (133) 

Then,  for  the  dimensional  mobility  K^m  and  diffusion  coefficient  Df]tm ,  we  hnd  the 
Einstein  relation 

kT 

Ddim  K dim  •  (134) 

e 

Moreover,  in  a  course  of  a  whole  non-stationary  process,  a  relationship  (130) 
between  ips  and  <pe  holds,  i.e.,  when  obtaining  the  solution  of  the  mobility  problem 
at  small  e  we  obtain,  also,  the  solution  of  the  diffusion  problem  at  small  5.  In  Report 
[28],  the  non-stationary  mobility  problem  was  solved  via  the  moment  method  and 
(pe(c,t)  was  built  up  for  several  interaction  models. 

Earlier,  the  stationary  solution  of  the  diffusion  problem  (fs(c,  oo)  for  HS-model, 
was  built  up  when  solving  an  integral  equation  in  [38] .  The  calculations  were  carried 
out  up  to  c  =  6.  If,  in  this  region,  using  a  formula  (130),  (ps(c,  oo)  was  expressed 
via  the  results  of  our  calculations  and  involve  the  corresponding  normalization,  they 
coincided  totally  with  the  results  [38]. 

Once  more,  it  is  confirmation  of  a  high  accuracy  those  data  which  we  obtain 
when  solving  the  Boltzmann  equation  via  the  moment  method  with  the  very  large 
number  of  the  MEs. 

The  study  of  diffusion  at  the  electric  held  was  performed  using  the  perturbation 
theory  via  the  two-stage  non-stationary  method. 

The  report  proves  a  possibility  to  use  the  perturbation  theory  for  the  set  of  the 
problem  being  under  consideration.  It  was  made  on  an  example  when  the  electric 
held  switch-on  developed  two  stage:  hrst,  a  sufficiently  strong  electric  held  e0  was 
switched  on,  second,  along  the  same  direction,  a  small  additional  electric  held  ej 
was  switched  on  being  considered  as  a  perturbation.  Such  a  test  turns  out  to  be 
true  as  there  is  another  mode  to  obtain  the  hnal  stationary  state  when  switching  on 
initially,  at  once,  a  held  eo  +  €\. 


At  the  switch-on  of  the  electric  held  e  along  z-axis,  the  moment  system,  describ¬ 
ing  the  non-stationary  process,  takes  a  form 


dCri  /  2  r(l  +  1) 

~df  +  e  V  2/ +  3 


(135) 
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When  the  system  being  solved,  we  find  the  dependences  Cri(t).  As  a  stationary 
solution  Cri( oo),  it  is  chosen  that  at  sufficiently  large  time  t.  With  these  moments, 
the  stationary  DF  is  restored  at  the  end  of  the  first  stage  of  the  process;  it  is  denoted 
by  /o(c). 

At  the  second  stage,  the  non-stationary  process  is  under  consideration  after  the 
switch-on  of  a  concentration  S  at  t  —  0.  The  DF  is  to  be  find  at  this  stage  as  follows 

f(c,t)  =  /0(c)  +  6M(ps(c,t),  (136) 


where  (ps(c,  0)  =  0. 

In  the  case  of  the  longitudinal  diffusion,  a  function  (pg(c,t)  is  axially  symmetric 
and,  in  the  moment  method,  it  is  expanded  in  terms  of  the  Legendre  polynomials. 
In  the  case  of  the  transverse  diffusion,  the  DF  losses  its  axial  symmetry  and  its 
expansion  is  carried  out  in  terms  of  the  spherical  harmonics  with  all  indices  l  and 
m.  The  moment  of  a  function  M(v)tps(c)  are  denoted  by  cri  and  cr;m,  respectively. 

After  rather  cumbersome  transformations,  the  moment  equations  for  second 
stage  were  deduced.  In  the  case  of  the  longitudinal  diffusion,  we  have 

~dT  +  "(27+3  r(/  +  ~  2tZ T^-1)  =  Y,Ar,n,lCrul- 


l  +  l 

2/ +  3 


21 

and,  in  the  case  of  the  transverse  diffusion, 


{(l  +  r  +  + 

_ ( /Wo)  / — *(o)  1  \ 

_  Y  1  W+lZ-l )  J  ) 


(137) 


da 


rim 


dt 


+  6 


2/ +  3 


r{l  +  l)cr_i;+im  — 


21-1 


lcrl— 


rl—lm 


r  1 


2/  +  3  (^  ^  r  ^  3/2)Wdi  +  i4-Wi)^i 

“srhOAh  -  c»H)^  (las) 

Here,  upper  index  0  of  the  moments  Cr[  denotes  that  those  are  the  moments  of  a 
function  /0,  i.e.,  the  DF  which  is  developed  at  the  end  of  the  first  stage,  being  the 
initial  state  for  the  second  one. 

When  increasing  the  electric  field,  a  function  ^(c),  in  the  case  of  the  longitudinal 
diffusion,  takes  the  very  large  value  in  a  domain  of  large  cz  and,  for  the  solution  to 
be  clear,  it  is  more  suitable  to  use  a  function  ^(c)  which  is  determined  from 


f{c,t)  =  /o(c)(l  +  S'lpsic,  t))- 


and  can  be  easily  expressed  via  935(c). 


(139) 
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a) 


b) 


Figure  18:  The  dependence  of  the  diffusion  tensor  on 
the  electric  field  e  for  the  various  interaction  models. 

Curve  1  -HS-model,  2  -  CEM-model,  3  -  CEHS- 
modcl.  Fig. a  -  longitudinal  diffusion  coefficient  D l. 

Fig.b  -  transverse  diffusion  coefficient  Dt- 

The  calculations  of  the  DF  evolution  were  carried  out  at  various  values  of  the 
electric  field  eo  and  three  interaction  models:  CEM-model,  HS-model,  and  CEHS- 


53 


model.  The  calculation  results  were  presented  in  details  in  the  11th  quarterly  re¬ 
port  [37]  in  which  there  are  many  figures  to  illustrate  the  evolution  of  both  <ps{c) 
(or  -05 (c))  and  the  coefficients  of  the  longitudinal  and  transverse  diffusion.  Here, 
we  restrict  ourselves  by  unique  figure  -  Fig.  18.  In  this  figure,  one  can  see  that 
the  stationary  coefficients  of  the  longitudinal  Dl  and  transverse  Dt  diffusion  are 
strongly  influenced  by  the  interaction  model. 

With  an  increase  in  the  electric  field,  for  all  interaction  models,  a  distinction 
between  Dl  and  Dt  increases.  The  more  pronounced  case  is  CEM-model.  Here, 
Dt  does  not  depend  on  e  at  all,  and  an  addition  in  Dl  increases  proportionally  to 
e2. 

For  HS-model,  the  transverse  diffusion  coefficient  increases  with  an  increase  in 
e.  However,  this  increase  lags  sufficiently  behind  an  increase  in  Dl  with  an  increase 
in  e.  In  the  case  of  CEHS-model,  Dt  decreases  with  an  increase  in  e. 

For  these  two  models  (HS-  and  CEHS-models),  the  dependences  T)l(c)  are  in¬ 
distinguishable  in  practice. 

The  general  results  obtained  during  Quarter  11: 

1.  To  calculate  the  diffusion  at  the  electric  field,  the  perturbation  theory  is 
developed.  Along  with  the  stationary  state,  the  non-stationary  process  developed 
after  the  switch-on  of  a  density  gradient  is  calculated. 

2.  Besides  the  diffusion  coefficient,  the  perturbation  DF  is  considered,  ft  is 
shown  that,  to  solve  the  problem  at  not  so  strong  electric  field,  the  standard  moment 
method  can  be  successfully  used  in  which  the  matrix  elements  of  large  indices  are 
involved. 

3.  Any  use  of  the  moment  method  and  perturbation  theory  is  confirmed  by  an 
example  when,  instead  of  the  gradient  there  is  a  small  extra  electric  field,  i.e. ,  under 
consideration  is  a  two  stage  process  of  the  electric  field  switch-on. 

4.  It  is  shown  that,  for  an  accurate  determination  of  the  perturbation  DF, 
the  high  accuracy  knowledge  on  the  stationary  DF  is  required  at  the  electric  field. 
As  it  is  shown  in  the  previous  reports,  at  the  strong  electric  field,  any  solution  of 
the  problem  under  consideration  via  the  standard  moment  method  meets  the  well 
defined  difficulties  related  with  the  Grad  criterion  violation. 

5.  When  obtaining  the  exact  stationary  solution  at  the  electric  field,  the  fine 
capabilities  of  the  developed  perturbation  theory  are  demonstrated. 

6.  The  calculations  of  the  DF  formation  after  a  switch-on  of  a  concentration 
gradient  for  several  ion-atom  interaction  models  are  performed.  For  the  first  time, 
the  DF  corresponding  to  a  perturbation  due  to  ion  diffusion  under  the  electric  field 
is  built  up. 

7.  The  calculations  of  the  coefficients  of  the  longitudinal  and  transverse  diffusion 
are  carried  out  at  the  electric  field. 
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12  Conclusion 


This  project  is  based  on  earlier  results  concerning  the  calculations  of  the  matrix 
elements  (ME)  of  the  collision  integral  of  the  Boltzmann  equation.  It  was  main¬ 
tained  that  there  is  a  set  of  relations  between  the  non-linear  moments  which  are  the 
recurrence  relationships.  A  code  package  was  developed  with  which  aid  both  the 
linear  and  non-linear  MEs  of  the  large  indices  can  be  found.  Earlier,  the  lack  in 
such  MEs  does  restrain  any  development  of  the  moment  method. 

To  solve  the  problem  stated  in  this  project,  the  non-stationary  moment  method 
is  advanced  in  which,  via  solving  the  non-stationary  system  of  the  moment  equation, 
the  evolution  process  is  studied  after  instant  switching-off  the  electric  field  of  taken 
value.  Further,  this  method  can  be  used  to  study  the  ion  behaviour  at  the  non- 
stationary  electric  field  of  any  form. 

The  left  sides  of  the  moment  equations  corresponding  to  the  differential  part 
of  the  Boltzmann  equation  are  introduced.  The  moment  equation  are  deduced  in 
the  global  basis.  Such  approach  is  called  the  standard  moment  method.  Also, 
the  equations  in  the  local  basis  are  built  up  when  the  parameters  of  the  weighted 
Maxwellian  depend  on  a  coordinate  and  time.  The  corresponding  method  is  called 
the  modified  moment  method. 

Considerable  attention  has  been  given  to  the  transfer  to  the  dimensionless  equa¬ 
tions.  The  standard  normalization  is  determined  in  which,  for  any  cross-  section  of 
ion-atom  interaction,  the  time  unit  is  unambiguously  determined.  Often,  the  non¬ 
standard  normalization  is  used  for  which  the  ion  mobility  equals  unity  in  the  limit 
of  infinitesimally  small  electric  field. 

The  codes,  with  which  aid  the  non-stationary  moment  method  is  realized,  have 
been  developed.  These  codes  are  for  all  problems  under  consideration  in  this  project. 
They  involve  the  reading  of  the  ME  matrix  from  the  previously  calculated  Table 
stored  in  a  special  file  for  taken  interaction  model.  Then,  the  evolution  of  the  DF 
moments  is  calculated  and,  using  them,  the  DF  itself  is  determined. 

The  interaction  models  is  advanced  for  which  the  analytical  solution  can  be 
built  up.  This  is  CEM-model  corresponding  to  the  resonance  charge  exchange  with 
a  cross-section  being  inversely  proportional  to  the  relative  velocity.  At  the  electric 
field  switch-on  at  no  magnetic  field,  the  analytical  solutions  are  built  up  both  for 
the  moments  and  the  DF. 

Using  the  analytical  solution  method,  the  detailed  study  of  the  standard  moment 
method  convergence  was  carried  out.  It  was  shown  that  using  the  non-stationary 
approach  is  capable  to  advance  significantly  further  in  building  up  the  stationary 
regime  as  compared  with  using  the  iteration  method  of  the  solution  of  the  stationary 
system  of  the  moment  equation. 

For  several  interaction  models,  the  detailed  calculations  of  the  DF  evolution  and 
the  physical  moments  are  accomplished  after  the  electric  field  switch-on. 
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It  was  carried  out  the  study  of  the  DF  and  physical  moments  in  crossed  mutually 
perpendicular  electric  and  magnetic  fields,  It  was  found  that  the  transient  non¬ 
stationary  process  is  accompanied  with  oscillations  of  a  current  along  the  electric 
held  as  well  as  the  Hall  current.  The  calculations  were  carried  out  for  the  Coulomb 
interaction.  At  the  weak  magnetic  held,  the  effect  of  ion  running-off  is  observed, 
being  rapidly  suppressed  with  an  increase  in  the  magnetic  held. 

The  non-stationary  moment  method  is  generalized  to  describe  the  ion  diffusion 
at  the  electric  held.  Two-stage  process  is  considered.  The  hrst  includes  the  non- 
stationary  process  after  the  electric  hied  switch-on.  The  second  is  involved  after 
yielding  the  stationary  state  at  the  beginning  of  which  a  small  concentration  gradient 
is  switched  on.  For  several  interaction  models,  the  coefficients  of  the  longitudinal 
and  transverse  diffusion  are  calculated. 

In  all  problems  under  consideration  in  this  project,  besides  the  transport  coeffi¬ 
cients  and  physical  moments,  the  DF  is  built  up.  For  the  hrst  time,  under  investi¬ 
gation  is  the  DF  evolution  after  switching  on  the  electric  held. 

The  studies  performed  result  in  a  clear  determination  of  the  limits  of  the  standard 
moment  method.  It  is  shown  that  a  poor  convergence  is  due  to  violation  of  the  Grad 
criterion.  The  paths  are  drafted  for  future  development  of  the  method  to  solve 
the  Boltzmann  equation.  Certain  success  can  be  reached  when  transferring  to  the 
modified  moment  method.  More  radical  solution  of  the  problems  is  the  expansion 
in  terms  of  the  Sonine  polynomials.  It  is  shown  that  the  realization  of  this  new 
method  is  possible  due  to  arising  the  new  method  of  calculation  of  the  MEs  of  the 
collision  integral. 

Both  the  standard  moment  method  and  the  mew  method  cited  above  can  be 
recommended  in  solving  a  whole  set  of  the  problems  of  the  gas  kinetic  theory  of 
gases  and  plasma.  Those  include  the  structure  of  shock  wave  as  well  as  the  study 
of  the  near-wall  layers  concerning  the  transport  problems  and  physical-chemical 
processes. 
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